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FAILURE  ANALYSIS  BY  STATISTICAL  TECHNIQUES  (FAST) 

1.0  INTRODUCTION 

FAST  is  a tool  developed  over  the  past  decade  because  of  the  need 
for  evaluating  the  nuclear  survivability  of  strategic  weapons  systems. 
These  weapons  systems  happened  to  employ  many  sites  built  almost 
identically.  This  led  to  early  recognition  that  a statistical  approach 
would  be  useful  - to  take  into  account  such  random  variations  as  those 
in  soil  properties  and  construction  quality. 

The  FAST  methodology  has  been  applied  mostly  to  in-place 
weapons  systems  for  the  purpose  of  evaluating  inherent  system  hardness 
or  the  benefit  of  hardness  improvements.  The  FAST  technique  has 
also  been  used  in  safety  analysis  of  ships,  and  is  applicable  to  earth- 
quake design.  In  general,  it  has  wide  potential  applicability  to  the 
hardness /survivability  evaluation  of  any  military  or  civilian  system. 

The  purpose  of  this  document  is  to  transmit  the  FAST  capability  to  other 
potential  users  so  they  can  employ  it  in  their  applications. 

This  report  contains  the  following  information  needed  to  enable  a 
prospective  user  to  understand  and  apply  the  FAST  techniques: 

• A description  of  the  analysis  formulation  and  of  the 
inputs  needed  to  perform  system  assessment  studies 

• Documentation  of  the  computer  code  and  its  operation 

• Demonstration  of  the  code  operation  through  documentation 
of  a sample  problem  which  has  been  designed  to  address 
the  key  aspects  of  the  methodology  and  its  practical 
application. 


1 . 

2.  0 FAST  METHODOLOGY 

This  section  describes  the  methodology  of  the  FAST  technique, 
the  analysis  formulation  and  system  modeling  considerations. 

2.  1 FAST  CONCEPT 

The  hardness  evaluation  or  assessment  of  a complex  system 
subjected  to  a hostile  environment  requires  calculation  of  the  probability 
of  response  of  each  component  to  the  hostile  environment,  determination 
of  the  probability  of  failure  of  each  component  for  that  response,  and 
the  combination  of  the  component  probabilities  to  obtain  the  failure 
probability  of  the  system.  The  FAST  code  has  been  designed  to  perform 
this  evaluation  in  a manner  such  that  parametric  sensitivity,  trade-off 
and  optimization  studies  can  be  readily  accomplished. 

Typically,  all  system  components  are  identified  and  catalogued. 

In  order  to  more  conveniently  catalogue  the  system,  its  components  can 
be  grouped  into  subsystems.  Then  the  possibility  of  'ach  component 
or  subsystem  directly  or  indirectly  contributing  to  system  failure  is 
ascertained.  The  individual  failure  mechanisms  and  fragilities  of  the 
components  are  determined  and  related  to  parameters  of  the  hostile 
environments. 

A diagram  describing  the  FAST  methodology  is  shown  in 
Figure  2-1.  The  figure  illustrates  how,  for  a simplified  system,  selected 
parameters  of  the  total  hostile  environment  may  be  critical  to  the  sur- 
vivability of  a hardened  system.  The  free  field  environments  are 
transformed  by  transfer  functions  to  establish  local  system  responses  to 
the  environments,  which  are  in  turn  used  to  predict  component  failure 
probabilities. 

In  the  FAST  code,  each  component  failure  probability  is  modeled 
by  a fragility  curve  that  defines  for  each  component  or  subsystem  the 
probability  of  failure  as  a function  of  the  local  system  response  to  the 
free  field  environments,  Component  probabilities  of  failure  are 
combined  in  system  network  equations  to  compute  system  probability 
of  survivability.  The  system  network  is  a functional  description 
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which  specifies  the  series /parallel  relationship  between  components. 
Components  are  in  series  if  all  components  are  required  to  accomplish 
an  essential  system  function,  ar.d  are  in  parallel  if  any  one  could  perform 
the  essential  function. 


An  important  facet  of  the  FAST  methodology  is  treatment  of  the 
underlying  uncertainties  in  predicting  system  hardness.  Of  the  four 
system/environment  inputs  to  FAST,  the  system  network  is  the  only 
one  that  must  be  known  exactly.  The  code  acknowledges  and  accommo- 
dates finite  uncertainty  in  modeling  environments,  transfer  functions  and 
component  fragilities.  Uncertainty  in  environment  estimates  is  often 
due  to  the  lack  of  adequate  analytic  or  empirical  models  for  scaling 
nuclear  weapon  effects.  This  is  also  true  of  the  uncertainties  ascribed 
to  transfer  functions.  Uncertainties  in  component  fragilities  are 
primarily  a consequence  of  insufficient  test  data  on  the  components  at 
levels  near  and  beyond  failure.  One  of  the  most  valuable  features  of 
FAST  is  that  these  uncertainties  are  individually  modeled  and  are 
properly  accounted  for  in  the  calculation  of  system  failure  probability. 
Two  categories  of  variations  are  recognized  by  the  FAST  code,  namely, 
random  and  systematic  variations. 

The  fundamental  difference  between  random  and  systematic  varia- 
tions is  that  a systematic  variation  extends  uniformly  over  a population 
of  facilities  whereas  random  variations  extend  non-uniformly  from  one 
facility  to  another.  Random  variations  tend  to  average  out  over  a large 
population  whereas  systematic  variations  do  not.  It  is  noted  that 
systematic  variations  are  usually  reducible  by  test  and/or  analysis 
programs  which  improve  models  of  the  phenomenology  and  system 
behavior.  Random  variations  may  also  be  reduced  under  some  conditions 
but  the  value  of  such  improvements  is  questionable  for  systems  with 
large  populations.  The  FAST  code  treats  both  random  and  systematic 
variations  of  environments,  fragilities,  and  transfer  functions.  The 
specific  treatment  of  variation  and  correlation  of  FAST  parameters 
is  noted  in  Table  2-A, 
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Table  2-A.  Variations  and  Correlations  of  FAST  Parameters 


Random 


Systematic 


Parameter  Variat 

Environments  Yes 

Transfer  Functions  Yes 
Fragilities  Yes 


Variation 

Correlation 

Variation 

Correlat 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

No 

Yes 

No 

Yes 

No 

Another  important  feature  of  FAST  is  that  the  code  accounts  for 
correlation  or  covariance  between  variables.  For  example,  if  two 
components  are  affected  by  strongly  correlated  environments,  they  will 
either  both  tend  to  fail  or  survive  with  the  probability  of  both  surviving 
being  higher  than  if  the  environments  are  uncorrelated.  Alternatively, 
if  the  environments  are  anti-correlated,  the  probability  of  survival  of 
the  system  can  be  substantially  less  than  for  the  uncorrelated  case. 
Relatively  strong  correlation  or  anti- correlation  between  environments 
is  not  unusual.  For  example,  a specific  system  might  be  vulnerable 
to  a combination  of  air  induced  ground  shock,  direct  induced  ground 
shock  and  debris.  Each  of  these  is  strongly  affected  by  soil  stiffness 
and  therefore  is  anti-correlated,  correlated  and  anti-correlated, 
respectively,  with  soil  stiffness.  FAST  treats  correlations  between 
both  systematic  and  random  variations  of  environment  parameters. 
Correlations  between  fragility  parameters  and  transfer  functions  are 
not  treated  except  for  azimuth  sensitivity  of  transfer  functions,  although 
such  correlation  is  easily  incorporated  into  the  code,  and  has  been 
incorporated  in  some  versions  of  FAST.  The  lack  of  use  of  these  options 
has  led  to  their  deletion  from  contemporary  versions  of  the  code. 

The  FAST  approach  combines  separate  estimates  of  the  environ- 
ments, transfer  functions,  fragilities,  corresponding  correlation  and 
systematic  variability  estimates,  and  a system  functional  network  to 
provide  probabilities  of  survival  of  the  system,  subsystems,  and 
components.  Probability  of  survival  is  given  as  a best  estimate  value 
along  with  certainty  (or  confidence)  bands  for  any  pre-selected  range 
from  a natural  or  manmade  energy  source. 
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The  FAST  code  is  standardized  for  systems  deployed  in  large 
numbers  and  built  almost  identically.  The  survivability  calculations 
are  performed  in  a fashion  to  discriminate  between  random  and  system- 
atic variations  in  the  survivability  statistics.  This  is  accomplished  by 
averaging  the  survivability  statistics  over  the  random  variations  to 
suppress  this  kind  of  variation  (but  account  for  non-linearities),  while 
preserving  and  displaying  the  variation  in  the  survivability  statistics 
due  to  systematic  variations  in  the  inputs.  However,  FAST  in  its 
current  form  can  also  consider  applications  involving  survivability  of 
systems  employing  only  one  or  a few  sites,  if  the  random  and  system- 
atic variations  can  be  combined. 


2.2  FAST  FORMULATION 


The  general  methodology  of  the  FAST  technique  allows  for  quite 
flexible  descriptions  of  parameters;  however,  the  existing  code  has 
very  specific  forms  of  the  environments,  transfer  functions,  fragilities, 
and  system  networks  which  have  been  found  to  be  the  most  useful  form 
for  code  application  by  the  user.  Specifically,  the  analysis  has  been 
developed  for  convenience  of  application  to  proliferated  systems  hardened 
against  nuclear  weapons  effects.  However,  parameters  can  be  given 
more  general  formulation  should  this  be  warranted  by  other  applications. 


The  computational  processes  built  into  the  current  code  are  not 
the  only  ones  that  could  be  used  for  the  FAST  technique;  however,  they 
have  allowed  for  the  most  general  employment  of  the  FAST  concept 
compared  with  other  computational  processes  used  in  the  past. 


2.2.1  Environments 


As  previously  discussed,  the  failure  probabilities  for  the  component 
fragilities  are  derived  from  free  field  hostile  environments  through 
transfer  functions.  The  FAST  hostile  environment  modeling  is  quite 
general.  For  example,  environments  may  be  variables  which  can 
directly  cause  failure  of  sensitive  components,  such  as  ground  shock, 
air  blast  or  debris  parameters.  Alternatively,  they  may  be  variables 
which  indirectly  contribute  to  failure,  such  as  depth  to  water  table  if 
failure  might  be  caused  by  flooding  after  some  structural  element  has 
failed. 


At  a reference  distance  from  a specific  nuclear  explosion,  a set 
ol  nominal  values  for  the  free  field  environments  is  defined.  These 
nominal  values  are  not  the  exact  values  that  would  occur  for  each 
explosion  if  a large  number  of  explosions  were  monitored.  The  varia- 
tions from  the  nominal  value  are  both  random  variations,  which  cover 
factors  like  topography  variations  at  the  reference  distance  from  a 
system,  or  geology  differences  from  one  site  to  another,  and  systematic 
variations,  which  deal  with  things  such  as  uncertainty  about  the  cratering 
efficiency  of  any  particular  bomb  in  a specific  geology  and  our  uncertainty 
in  defining  the  specific  geology  at  the  sites  of  interest.  A set  of  nominal 
values  may  be  calculated  by  predicting  a set  of  best  estimate  environ- 
ments for  each  site  in  the  group  (or  each  site/azimuth)  and  then  statis- 
tically calculating  the  mean  value  for  each  environment  in  the  set: 
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(2-1) 


This  set  of  nominal  values  is  referred  to  as  the  "vector  of  environment 
means"  and  is  identified  in  FAST  as  ULN.  These  are  of  course  related 
to  specific  reference  threat  conditions.  An  example  giving  median 
environments  for  a 600  psi  overpressure  pulse  from  a 5 MT  weapon 
follows. 

Each  environment  parameter  is  represented  statistically  to  take 
into  account  the  random  variations,  and  systematic  variations,  as 
noted  in  Section  2.  1.  The  normal  or  log-normal  distribution  is  used  to 
model  each  environment.  An  example  of  the  normal  distribution  is 
shown  in  the  following  sketch.  In  the  log-normal  distribution,  the 
logarithms  of  the  parameters  a^e  distributed  normally.  The  log-normal 
distribution  for  a single  variable  is  given  by 


4>  (log  x) 


n/2tt  x tr(log  x) 
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FAST  Code 


Parameter 

Designation 

Example 

Reference  Threat 

NOM  RWH 

Range 

R 

o 

3000  feet 

Yield 

w 

o 

5 MT 

Height  of  Burst 

H 

o 

3 feet 

Vector  of  Means 

ULN 

Air  Induced 

Displacement 

E01 

15  inches 

Velocity 

E02 

380  in. /sec 

Acceleration 

E03 

60  g's 

Direct  Induced 

Displacement 

E04 

27  inches 

Velocity 

E05 

29  in. / sec 

Debris 

E06 

12  inches 

Water  Table 

E07 

50  feet 

Note:  The  overpressure/range  relationship  is  from  Ref.  (3). 


where  x is  the  environment  parameter,  cr(log  x)  is  the  standard 
deviation  and  <b  is  the  probability  density.  An  example  of  the  log-normal 
distribution  is  shewn  in  the  following  sketch.  Frequently,  strictly 
positive  random  variables  follow  the  log-normal  distribution  which  has 
proved  to  be  adequate  for  all  past  calculations;  however,  the  distributions 
could  be  expanded  to  include  other  density  functions. 


MEAN  - 0 

STANDARD  DEVIATION 


MEAN  = I 
K-FACTOR  - Z 


NORMAL  DISTRIBUTION 


LOG-NORMAL  DISTRIBUTION 
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As  discussed  in  Section  2.  1,  the  variability  of  different  environments 
is  not  usually  independent,  so  that  correlations  between  environments 
must  be  properly  modeled  in  the  code.  Correlation  is  modeled  by 
correlation  coefficients  which  mathematically  are  defined  by 


P, 


Cov  (x,  y) 


O-  (T 

x y 


(2-3) 


where  cr  and  <r  are  the  standard  deviations  of  the  random  variables  x 
x y 

and  y,  and  Cov  (x,  y)  is  defined  as 


Cov  (x,  y)  = J J" (x-x)  (y— y ) f(x, y ) dxdy  (2-4) 


The  quantities  x and  y are  the  mean  values  of  x and  y,  and  f (x,  y)  is 


the  joint  probability  density  function.  When  px^  is  zero  the  random 


variables  are  linearly  independent;  when  px^  is  unity  the  random  vari- 
ables are  linearly  correlated;  and  where  p is  minus  one  the  random 
7 rxy 


variables  are  linearly  anti-correlated.  For  multivariate  distributions, 
the  variances  and  correlation  coefficients  can  be  expressed  in  matrix 
form  as  shown. 


2 


P12V2 


P,  O',  (T 
L In  1 n 


P1 2 °”l  °2 PlnVn 

S P2na2  °n 

cr'2 


(2-5) 


This  array  is  known  as  the  covariance  matrix.  The  o's  are  the  standard 
deviations  of  the  environments  (linear  or  log)  so  that  the  matrix  can 
be  used  to  define  the  joint  probability  density  function  for  all  the 
environments  as  shown. 


4,(xr  x2,  xn) 


(x  - x)2  (x  - x) 


(2-6) 


1.3 
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where  (x  - x)  is  an  n-component  column  vector  comprising  all  of  the 
environment  parameters,  (j T^lf)  is  its  transpose,  D"*  is  the  inverse  of 
the  covariance  matrix  and  | E | is  the  determinant  of  X.  It  is  noted  that 
some  of  the  variables  in  this  distribution  may  be  logarithms  of  environ- 
ment parameters,  while  others  are  the  environment  p.  'ameters 
themselves.  The  only  requirement  is  that  the  distribution  of  the 
variable,  or  its  logarithm,  be  normal.  References  1 and  2 contain 
more  detailed  information  on  the  properties  of  multivariate  and  log- 
normal distributions.  Two  covariance  matrices  are  modeled  for  each 
FAST  calculation.  One  corresponds  to  the  random  variation,  and  the 
other  to  the  systematic  variation.  The  former  is  input  as  the  S-matrix 
->nd  can  be  derived  from  best-estimate  calculations  for  each  siting 
area  (or  azimuth/site  group).  The  second  is  described  by  a K-factor 
vector,  defining  the  systematic  variation  in  each  environmental  estimate, 
and  an  A-matrix  characterizing  the  correlation  of  these  variations  (see 
Section  3). 

These  covariance  matrices  are  not  scaled;  however,  individual 
matrices  can  be  specified  for  each  threat  condition.  Experience 
indicates  that  the  matrices  for  weapon  environment  parameters  are  at 
least  approximately  constant  if  the  scaling  laws  are  linear  (in  their 
logarithmic  form). 

The  multivariate  normal  distribution  is  defined  by  its  mean 
vector,  x,  and  its  covariance  matrix,  E.  The  role  of  each  of  these  in 
FAST  analysis  will  now  be  discussed. 

Median  environment  parameters  are  first  defined  for  a reference 
threat  condition.  If  the  analysis  concerns  a facility  subjected  to  nuclear 
attick,  the  median  environment  parameters  are  defined  for  a baseline 
threat  condition  of  yield  (W^),  miss  distance  range  (RQ)  or  corresponding 
mean  overpressure  (Pq),  and  height  of  burst  (Hq).  These  environments 
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are  then  scaled  to  the  other  threat  condition  of  interest  (W,  R or  P 
and  H),  using  either  or  the  two  forms  of  scaling  laws  shown 


(2- 7a) 


or 


(2-  7b) 


where  a , (3  and  V are  scaling  factors,  and  Eq  and  E are  the  baseline 
and  scaled  environment  parameters,  respectively. 

The  FAST  code  contains  the  Brode  overpressure  relationship 
with  height  of  burst  and  range,  specified  in  Reference  3,  Through 
this  equivalence  between  nominal  overpressure  and  range,  the  choice 
of  scaling  law  for  each  environment  parameter  is  purely  a matter  of 
convenience. 

The  height  of  burst  scaling  is  designed  to  model  crater-related 
effects  for  near  surface  burst  conditions.  A much  more  complex  scaling 
would  be  needed  to  cover  all  conditions  from  surface  to  above  optimum 
height  of  burst.  High  altitude  burst  conditions  can  be  considered  by 
suppressing  the  height  of  burst  scaling  and  letting  the  baseline  threat 
contain  the  height  of  burst  of  interest. 

2.2.2  Transfer  Functions 

Transfer  functions  relate  the  free  field  environments  to  the  local 
system  response  of  components.  In  normal  usage  a transfer  function 
will  take  a single  environment  (for  example  overpressure)  and  convert 
it  into  a response  which  governs  a system  failure  mode  (for  example 
shear  force  in  a facility  roof).  The  transfer  function  can  be  used  to 
account  for  structural  dynamics,  impedance  mismatcn,  impulse  effects, 
etc.  Moreover,  transfer  functions  can  be  used  in  ways  w'hich  simplify, 
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clarify,  or  replace  environment  generation.  Indeed,  environments, 
transfer  functions  and  fragilities  are  so  uniquely  correlated  that  any 
particular  detail  of  the  modeling  might  be  incorporated  in  any  one  as 
well  as  the  other. 


Transfer  functions  can  be  a function  of  the  attack  azimuth  (0). 
Problems  where  applied  loads  vary  with  azimuth,  due  to  non-axisymmetric 
hardware  or  site  conditions,  can  be  modeled  in  this  manner.  In  its 
simplest  form,  a transfer  function  (TF)  relates  ar.  environment  parameter 
(F)  to  a response  (R),  viz.  ; 


R = E x TF  (0) 


(2-8) 


A transfer  function  may  be  less  than  unity,  representing  attenuation, 
or  it  may  be  greater,  representing  amplification. 

Transfer  functions  are  represented  statistically  in  the  FAST  code. 
Systematic  variation  is  assumed  to  have  a beta  distribution  specified  by 
a K-iactor,  with  95%  of  the  distribution  falling  between  (median/K)  and 
(K  x median)  and  2-1/2%  at  each  end.  The  code  selects  the  best  fit 
from  the  nine  choices  of  the  beta  distributions  shown  in  Figure  2-2. 

This  table  can  be  replaced  to  provide  a different  family  of  distributions 
if  such  is  appropriate  to  a given  application.  As  discussed  in  2.  2.  3, 
t ho  same  table  is  used  for  fragility  systematic  variation. 

Random  variation  of  transfer  functions  is  introduced  into  the  FAST 
code  to  account  for  azimuth  sensitivity  of  components  in  the  analysis  of 
facilities  subjected  to  nuelea"  attack.  This  is  accomplished  by  use  of 
azimuth  sensitivity  tables.  An  azimuth  is  selected  at  random  (uniform 
d ist  r i out  ion)  and  used  for  all  components. 


In  the  more  general  form  of  the  transfer  function,  a single  system 
response  can  be  written  as  a linear  combination  of  up  to  30  environment 


parameters,  viz.; 


R 


E , x TF.  (0)  f E,  x TF,  (0)  + 

11  t-  C. 


(2-0) 
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Such  a formulation  is  very  useful  for  modeling  dynamic  response  when 
modal  participation  is  a factor. 

2.2.3  Fragilities 

A fragility  defines  the  probability  of  failure  of  a component  as  a 
function  of  a local  system  response  {frequently  a critical  local  environ- 
ment). In  many  applications,  failure  is  assumed  to  occur  if  the  component 
experiences  damage  sufficient  to  prevent  (by  itself  or  in  combination  with 
other  failures)  successful  completion  of  the  system  mission.  No  other 
damage,  no  matter  how  severe,  should  be  considered  failure.  However, 
the  code  is  simply  a tool,  and  if  another  failure  definition  is  more 
appropriate  for  a specific  application,  it  could  of  course  be  used. 

Possible  forms  of  fragilities  are  shown  in  Figure  2-3.  The  ordinate 
is  probability  of  failure,  ranging  from  0 to  1,  and  the  abscissa  is  a 
component  response  (ordinarily  a critical  local  environment)  which  can 
cause  failure.  A component  may  have  more  than  one  fragility,  which 
might  correspond  to  different  failure  modes  or  other  system  modeling 
considerations.  Each  fragility  must  represent  the  best  engineering 
estimate  of  the  actual  failure  probability  function,  not  a conservative 
estimate.  It  must  represent  as  much  likelihood  of  being  too  low  as 
being  too  high.  All  fragilities  are  treated  as  being  uncorrelated. 

A fragility  curve  is  a cumulative  probability  function  which,  in 
general,  is  sloped  due  to  postulated  differences  among  a group  of  like 
components  or  other  random  variations.  For  example,  individual 
components  built  to  the  same  drawings  may  fail  at  different  response 
levels  because  of  unintentional  differences  in  tolerance,  material 
properties,  parts  production  variation,  welding,  manufacturing  pro- 
cedures, maintenance,  operating  environment,  or  other  factors. 
Furthermore,  fragility  testing  considers  only  an  approximation  of 
the  local  environment,  while  fragility  analyses  consider  only  an 
approximate  model  of  the  component.  The  resulting  uncertainties  lead 
to  both  random  and  systematic  variations. 
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Figure  2-3A  shows  a cookie  cutter  fragility,  i.e.  , without  random 
variation;  whereas  2-3B  shows  a fragility  with  random  variation. 

Analysis  and/or  specific  tests  may  be  used  to  determine  the  steepness 
and/or  straightness  of  the  curve.  FAST  allows  the  general  fragility 
curve  to  be  composed  of,  at  most,  five  linear  segments  beginning  at 
zero  and  ending  at  one.  The  format  is  illustrated  in  Figure  2- 3D. 

The  upper  and  lower  dashed  boundaries  in  Figures  2-3C  and2-3D 
reflect  the  systematic  variation  or  uncertainty  inherent  in  establishing 
the  best  estimate  curves.  The  boundary  lines,  in  general,  need  not  be 
parallel  or  symmetric  about  the  nominal  curve. 

The  fragility  systematic  variation  within  the  uncertainty  band  is 
assumed  to  have  a beta  distribution  with  the  best  estimate  centered  at 
the  mode  of  the  distribution.  The  code  selects  the  best  fit  from  the 
nine  choj~.es  of  the  beta  distribution  shown  in  Figure  2-2. 

Azimuth  sensitivity  of  fragilities  can  also  be  accounted  for  in 
FAST.  The  simplest  example  of  azimuthal  effects  on  fragilities  is  a 
shock-isolated  component  with  different  amounts  of  rattlespace  in 
different  directions.  One  assigns  a reference  direction  for  the  problem 
and  all  azimuth  sensitivities  are  measured  from  the  reference  azimuth. 
Then  during  calculation,  random  azimuths  are  selected,  and  for  each 
azimuth  a’l  environment  parameters  are  converter  to  fragility  responses, 
takinginto  account  azimuth  sensitivity,  and  compared  with  critical  fragil- 
ity levels.  In  this  way  systems  which  have  the  same  weak  directions  for 
multiple  components  will  not  have  multiple  kills  traceable  to  weak 
direction  effects,  when  all  the  azimuthal  failures  occur  either  concur- 
rently or  not  at  all. 

2.2.4  System  Network 


The  logic  enabling  the  code  to  combine  component  failure 
probabilities  into  subsystem  and  system  failure  probabilities  solves 
probability  equations  based  upon  the  system  network.  The  code  uses 
two  levels  of  networks:  components  are  combined  to  form  subsystems 
and  subsystems  and/or  components  are  combined  to  form  systems. 
Figure  2-4  shows  several  examples  of  such  networks:  2-4A  is  a 
subsystem  of  three  components  in  series,  2-4B  a subsystem  of  three 
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components  in  parallel  and  2-4C  is  a subs/stem  including  both  series 
and  parallel  components.  Systems  are  defined  as  combinations  of 
subsystems  and/or  components  in  an  analogous  way,  as  indicated  by 
Figure  2-4D. 

For  components  in  series,  a failure  of  any  one  of  them  will  cause 
subsystem  failure.  Thus,  if  the  probabilities  of  failure  for  the  components 
in  Subsystem  1 of  Figure  2-4  are  P^qj»  Pjq2  anc*  ^103  *^en  the  proba- 
bility that  this  subsystem  fails  is 


SSI 


1 -<UP101><‘-P102)(1-P103> 


(2- 10a) 


For  components  in  parallel,  a failure  of  all  is  required  for  subsystem 
failure.  Thus,  the  probability  of  failure  for  Subsystem  2 is 


P = P . p • p 

*SS2  * 201  * 202  * 203 


(2-  10b) 


Any  subsystem  or  system  can  be  modeled  as  combinations  of  series 
and  parallel  elements.  Consequently,  the  failure  probabilities  for  the 
other  networks  of  Figure  2-4  are  given  by 


PSS3  = 1_  (1_P302  P303)(1'P301)] 


ri-(1-P304>  <1-p305>]  ><1**P306> 


(2-10c) 


PSYS1  = 1_  (1_PSS3  PSS4)  (i'PSSl)  (1‘PSS2)  (2-10d) 
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Derivation  and  programming  of  such  mathematical  expressions 
is  not  practical  for  complex  systems,  so  the  FAST  code  evolution  has 
developed  a simplified  automated  technique  for  this  purpose.  The  tech- 
nique employs  software  called  the  Boolean  Network  Compiler.  The 
relationship  between  the  various  components  and  subsystems  is  input, 
and  the  compiler  sets  up  the  correct  mathematical  operations  to  com- 
pute subsystem  and  system  survivability.  To  illustrate  the  simplicity 
of  this  technique,  equations  2 -10a  through  2-10d  would  be  input  as 
follows: 


= C101  + C102  + C103 


C201  * C202  * C203 


(2 -10a') 
(2-1  Ob’) 


= (C302  * C303  + C301)  * (C304  + C305)  + C306  (2-10c ') 


SYS1  = SSI  + SS2  + SS3  * SS4 


(2-1 Od ' ) 


The  Boolean  operators  + and  * denote  the  series  and  parallel  relation- 
ships, respectively. 

The  ability  to  name  several  systems  in  a single  run  by  varying 
one  or  more  subsystems  or  components  is  called  "ganging",  and  is 
useful  for  parametric  studies.  This  is  because  once  the  component 
failure  probabilities  have  been  computed,  they  can  be  combined  in 
different  ways  according  to  the  network  equations  to  obtain  various 
system  failure  probabilities, 

2.  2.  5 Survivability  Statistics 

This  section  describes  how  the  FAST  code  combines  the  statis- 
tical representations  of  environments,  fragilities,  and  transfer  func- 
tions in  accordance  with  the  system  network  definitions  to  obtain  a 
statistical  description  of  the  system  survivability. 
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2.  2.  5.  1 Survivability  Calculation 


The  FAST  code  is  standardized  for  proliferated  hardened  systems 
as  discussed  in  Section  2.  1.  Random  variables  affecting  system  surviv- 
ability are  appropriately  averaged  by  the  process,  while  the  system- 
atic variations  are  propagated  through  the  calculations  to  indicate  the 
level  of  certainty  (or  confidence)  in  the  system  survivability  result. 

The  Monte  Carlo  technique  is  used  in  performing  the  survivability 
calculations.  This  technique  samples  from  the  random  and  systematic 
distributions  of  the  environments,  fragilities,  and  transfer  functions 
to  obtain  values  used  to  compute  a system  survival  probability.  This 
process  is  repeated  many  times  to  obtain  a sample  of  survival  proba- 
bilities which  can  be  summarized  and  displayed  as  system  survivability 
statistics. 

An  overview  of  the  computation  process  is  diagrammed  in  Fig- 
ure 2-5  for  a simplified  system  of  just  two  -'omponents.  These  com- 
ponents, together  with  their  associated  environments,  fragilities,  and 
transfer  functions,  are  denoted  by  the  subscripts  1 and  2.  The  FAST 
code  accomplishes  the  sampling  process  in  two  stages,  to  discriminate 
between  random  and  systematic  variations. 

The  first  stage,  accomplished  in  an  outer  loop  of  the  code,  selects 
bias  values  from  the  systematic  variation  distributions  for  the  environ- 
ments, fragilities,  and  transfer  functions.  This  is  illustrated  by  the 
point  © (ej,  e2),  selected  from  the  environment  systematic  variation 
distribution  in  the  upper  left  in  Figure  2-5,  by  the  selection  of  the 
solid  lines  labeled  Biased  TFj  and  Biased  TF£  from  the  transfer  func- 
tion systematic  variation  distributions  and  by  selection  of  the  solid  lines 
labeled  Fragility  1 and  Fragility  2 from  the  fragility  systematic  variation 
distributions. 

The  systematic  variation  information  defined  in  the  outer  loop  is 
passed  to  the  inner  loop  where  the  second  stage  of  the  sampling  process 
is  accomplished.  Values  are  selected  from  the  random  variation  dis- 
tributions for  the  environments  and  transfer  functions  (azimuth  sensi- 
tivity). The  environment  values  are  chosen  from  the  multivariate 
normal  distribution  defined  by  the  mean  ((ej,  e^)  in  the  illustration)  and 
the  random  covariance  matrix  (Equation  2-5).  One  such  selection  is 
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indicated  by  the  symbol  * in  the  inner  loop  environment  distribution 
diagram  of  Figure  2-5.  If  the  environment  parameters  are  log- 
normally  distributed,  antilogs  are  taken  at  this  point. 

If  azimuth  sensitivity  tables  have  been  provided  for  any  of  the 
transfer  functions,  a single  random  number  is  used  for  selecting  from 
all  of  the  azimuth  tables.  Randomization  of  the  transfer  function  for 
Component  2 is  indicated  in  Figure  2-5. 

When  randomization  on  the  environment  parameters  and  transfer 
functions  is  complete,  calculation  of  component  probabilities  of  failure 
is  accomplished  as  indicated  by  the  nomographs  of  Figure  2-5.  The 
fragility  response  is  compufed  as  follows: 


Ri  = TFj  exp  (ej) 


(2-11) 


The  response  parameter  , R^,  is  input  to  the  fragility  curve  to  obtain 
the  component  probability  of  failure. 


Next  the  component  probabilities  of  failure  are  combined,  using 
the  logic  network  compiler  to  obtain  subsystem  and  system  probabilities 
of  failure. 


The  inner  loop  process  is  repeated  until  sample  size  or  conver- 
gence criteria  (discussed  in  the  following  section)  arc  satisfied.  The 
mean  probability  of  failure  is  computed  from  the  accumulated  inner 
loop  sample  statistics  for  every  component,  subsystem  and  system. 
These  sample  means  constitute  the  output  from  a single  outer  loop 
iteration.  The  process  is  repeated  for  many  outer  loop  iterations, 
accumulating  the  probability  of  failure  data  (inner  loop  sample  means) 
in  histogram  format.  This  process  continues  until  outer  loop  sample 
size  criteria  (discussed  below)  are  satisfied.  The  output  from  the 
outer  loop  iterations  forms  the  survivability  statistics. 
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2.  2,  5,  2 Convergence  Criteria 

The  data  generated  by  the  inner  and  outer  loops  can  be  represented 
by  an  array  as  follows. 


Row 

Averages 


P11  P12  ‘ ‘ Plm. 


P2 1 P22  ‘ * P2m. 


pl 

Po 


(2-12) 
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where  the  pij  are  failure  probabilities.  Each  row  in  this  array  contains 
the  probability  of  failure  data  for  a system  obtained  from  repeated 
iterations  of  the  inner  loop  for  one  iteration  of  the  outer  loop.  Each 
row  is  averaged  to  suppress  the  random  variation  within  the  row, 
yielding  the  p vector  shown  at  the  right  above,  whose  elements  are  a 
sample  from  the  distribution  of  system  probability  of  failure  perturbed 
by  bias  variations  on  the  inputs. 


The  FAST  code  contains  an  algorithm  to  control  the  sampling 
accuracy  from  the  random  and  bias  distributions.  The  objective  is  to 
minimize  computer  cost  without  sacrificing  necessary  accuracy. 

If  too  few  iterations  are  done  on  the  inner  loop,  the  variance  of 
the  average  probability  of  failure  will  be  excessive  and  appear  in  the 
bias  distribution  of  the  outer  loop,  which  is  an  undesirable  condition. 
To  control  this,  data  on  variability  of  system  probability  of  failure 
within  the  inner  loop  is  collected  over  the  first  5 (J  outer  loop  iterations 
to  establish  the  number  of  inner  loop  iterations  required  to  give  the 
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specified  accuracy.  The  inner  Loop  convergence  algorithm  makes  use 
of  the  following  equation. 
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(2-13) 
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£i  L (Ni  -1* 


where 


Nx(k) 


Estimated  number  of  inner  loop  iterations  required 
after  k (150)  outer  loop  iterations  (Note:  NTT(k)-5 
and  Nt  (k)  s 2Nj  (k-1)  are  required). 


N; 


Number  of  inner  loop  iterations  that  were  done  in 
the  Ik  outer  loop  iteration 


Number  of  outer  loop  iterations  performed  (k  s 50). 


Pij 


Probability  of  failure  from  the  jt*1  inner  loop  itera- 
tion of  the  i*h  outer  loop  iteration 


Maximum  error  bound  on  the  inner  loop  average 
probability  of  failure 


Equation  (2-13)  ? ^ applied  repeatedly  during  the  first  50  outer  loop 
iterations.  After  50  outer  loop  iterations,  the  number  of  inner  loop 


iterations  is  fixed  at  the  final  value  of  N.j,(50). 


Convergence  of  the  outer  loop  is  not  tested  until  fifty  iterations 
have  been  performed.  At  this  time,  the  following  computation  is 
performed, 
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where 

N(k)  = Estimated  required  number  outer  loop  iterations 
k = Number  of  outer  loop  iterations  performed 

2 

= Required  error  bound  on  the  variance  of  the  outer 
loop  average  probability  of  failure. 

= Mean  of  the  i^  row  in  Eq.  (2-12) 

If  N(k)  is  moie  than  the  number  of  outer  loop  iterations  that  have  been 
performed,  additional  iterations  are  done  in  increments  of  50  and  the 
check  is  repeated  until  N(k)  is  less  than  the  number  of  outer  loop 
iterations  performed. 

After  convergence  on  the  mean  is  achieved,  an  additional  test 
for  convergence  on  the  median  is  performed.  This  test  uses  the 
following 


where 


f(Ck  + \ kC(l-C))  - f(Ck  - \ kC(l-C))  < 2c02 


Number  of  outer  loops  performed 


(2-15) 


C 0.  5,  for  the  median 

f(x)  - Fraction  of  the  elements  in  the  p vector  of  array 
(2-12)  which  are  less  then  x (1  s x s k) 

- Required  accuracy  on  the  median  probability  of 
failure. 

If  the  inequality  is  satisfied,  outer  loop  iteration  stops. 

The  interpretation  of  (2-15)  is  as  follows: 

Given  that  k outer  loops  were  performed,  the  p vector  of  array 
(2-12)  is  arranged  in  order  of  increasing  magnitude.  The  resultant 
reordered  array  is  called  f(x),  where  x is  the  new  index  (low  index, 
low  magnitude). 
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Then  equation  (2-15)  says  that  the  average  quantitative  difference 
between  two  values  of  f,  one  on  each  side  of  the  median,  and  the  com- 
puted median  is  a measure  of  the  uncertainty  in  the  computed  median. 
Additional  outer  loops  are  performed,  in  increments  of  50,  until  the 
1 inequality  (2-15)  is  satisfied  or  until  a maximum  of  625  outer  loops 

have  been  performed. 

The  number  of  outer  loop  iterations  required  by  the  last  criterion 
can  become  large  as  the  average  of  the  outer  loop  failure  probability 
distribution  approaches  C = 0.  5.  This  is  the  reason  that  the  upper 
limit  of  625  outer  loop  iterations  has  been  imposed.  This  upper  limit 
can  be  changed  if  a different  value  is  appropriate  for  a given  application. 

2.  2.  5.  3 Environment  Parameter  Vector  Generation 

The  F4.ST  code  requires  large  numbers  of  vectors  of  environment 
parameters  drawn  from  the  distributions  of  systematic  and  random 
variations.  The  distributions  are  multivariate  normal  (after  taking 
logarithms  of  log-normal  environment  parameters)  with  mean  p and 
covariance  Z.  To  provide  computational  efficiency  and  simplicity  of 


programming,  both  the  systematic  and  random  covariance  matrices 
are  operated  upon  by  a matrix  decomposition  approach.  It  is  noted 
that  the  approach  used  in  FAST  requires  that  the  matrix  Z be  positive 
semi-definite  (most  routines  of  this  type  require  Z to  be  positive 
definite).  This  means  that  the  code  can  accommodate  analysis  formu- 
lations wherein  some  environments  are  linear  combinations  of  others. 
This  section  summarizes  the  theory  of  the  matrix  decomposition 
formulation, 

A matrix  Z (which  must  be  positive  semi-definite)  can  be  decom- 
posed into  a matrix  C such  that 

C C = Z (2-16) 
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where  C is  the  transpose  of  C.  Then  the  required  random  vector  x is 
calculated  from  a vector  y (whose  elements  are  independently  and 
normally  distributed  random  variables  with  zero  means  and  unit 
variances)  as  follows: 


x = p.  + Cy 


(2-17) 


To  see  that  x has  the  desired  distribution,  assume  CC  = I and 
let  x'  = Cy,  where  y is  a vector  of  uncorrelated  random  variables 
with  mean  zero  and  unit  variance  (and  hence,  E [yyj  = I,  where 
E {•}  denotes  expectation).  Then 


E x'}  = 


- EiCyy  C} 
= C E {y  y }C 

/v 

= c i r 

/V 

- C C 


(2-18) 


Thus,  x'  has  zero  mean  and  covariance  2.  Finally,  setting  x = p + x' 
completes  the  demonstration. 

The  method  for  computing  C from  2,  known  as  the  square  root 
method,  is  given  in  Reference  (4),  which  should  be  consulted  for  ie 
basic  theory.  Most  of  the  algorithms  in  the  literature  require  that  2 
be  positive  definite.  The  matrix  decomposition  algorithm  in  this 
analysis  incorporates  a modification  which  relaxes  that  requirement, 
so  that  2 canbe  positive  semi-definite  (i.  e,  , its  rank,  r,  can  be  less 
than  its  size,  n).  In  the  modified  routine,  each  successive  pivotal 
element,  C (k,k),  is  selected  as  the  largest  remaining  leading  diagonal 
clement  rather  than  in  numerical  order  as  in  the  standard  algorithm. 
Thus,  if  2 has  rank  r<n  (the  size  of  2)  after  r applications  of  the 
algorithm,  the  remaining  leading  diagonal  elements,  C (k,k),  will  be 
approximately  zero  (i.  e.  , equal  to  zero  except  for  computational 
errors  due  to  rounding).  Rows  containing  diagonal  elements, 

C (k,k)  0,  are  set  equal  to  zero.  After  the  decomposition  routine 

has  been  completed,  the  matrix  contains  r non-zero  diagonal  elements 
and  n-r  zero  diagonal  elements.  Then  rank  of  C (and  2)  is  r. 
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2.  2.  5.  4 Output  Displays 


Results  of  FAST  calculations  can  be  displayed  in  a variety  of 
fashions  as  illustrated  in  Figure  2-6.  One  important  display  is  that  of 
the  system  survivability  together  with  associated  uncertainty,  as  a 
function  of  overpressure  or  miss  distance.  (See  upper  left  of  Fig- 
ure 2-6).  In  the  upper  right  of  this  figure  is  plotted  the  median  value 
of  the  probability  of  survival  of  the  system  and  its  subsystems,  as 
functions  of  overpressure.  This  output  format  is  valuable  for  identi- 
fying subsystems  which  are  weak  link  items,  and  hence  potential 
candidates  for  hardening. 

The  other  three  diagrams  of  the  figure  each  show  the  median 
probability  of  survival  of  a critical  subsystem,  together  with  the 
medians  of  the  components,  all  as  functions  of  overpressure.  This 
format  enables  identification  of  components  which  are  weak  link  items 
in  the  subsystem.  Taken  together,  these  displays  enable  identification 
of  those  subsystems  and  components  contributing  most  to  the  failure 
of  the  system. 

2.  3 SYSTEM  MODELING  CONSIDERATIONS 

The  derivation  of  inputs  for  probabilistic  assessments  of  a 
system's  capability  to  function  after  being  subjected  to  hostile  environ- 
ments can  proceed  in  several  ways.  No  single  approach  is  best  for 
all  systems,  as  some  problems  naturally  emphasize  environments, 
other  problems  highlight  transfer  functions  or  fragilities.  It  is 
necessary  to  take  a broad  view  prior  to  developing  the  mathematical 
model  of  the  system,  since  the  approach  which  is  taken  in  formulating 
the  problem  will  influence  the  system  network. 

The  first  consideration  is  to  specify  the  definition  for  system 
failure  or,  conversely,  system  survival.  A system  is  generally  said 
to  survive  if  it  retains  the  capability  to  perform  its  mission  for  a 
stated  period  of  time  after  being  subjected  to  one  or  more  applications 
of  hostile  environments.  Some  systems  may  have  multiple  missions, 
such  as  both  tactical  and  strategic,  each  of  which  utilizes  different 
combinations  of  subsystems.  Another  definition  of  system  survival 
(repairable  survival)  might  allow  a system  to  be  inoperable  for  some 
stated  time  interval  provided  that  it  could  be  returned  to  service 


33 


(\  e.  , regain  the  capability  to  perform  its  mission)  within  the  time 
interval.  Obviously,  the  differences  between  single  and  multiple 
attack  conditions  and  emergency,  extended  and  repairable  survival 
definitions  will  greatly  affect  the  derivation  of  inputs  for  FAST  code 
system  simulations. 

This  section  considers  various  techniques  and  practical  considera- 
tions for  developing  the  system  network,  fragilities,  environments,  and 
transfer  functions.  To  do  this,  the  following  data  needs  to  be  considered: 

a.  Environments 

• Median  values  for  baseline  threats. 

• Random  and  systematic  variation  covariance  matrices 
for  baseline  threats. 

• Scaling  parameters 

b.  Transfer  functions  relating  environments  with  corre- 
sponding fragility  parameters 

• Best  estimate  values 

• Uncertainty  bounds 

c.  Fragility  curves  for  each  network  element 

• Best  estimate  curve 

• Bounding  curves 

d.  System  network. 


2.  3.  1 Environments 

As  previously  discussed,  the  FAST  code  accepts  a probabilistic 
definition  of  environments,  with  the  multivariate  log-normal  distribu- 
tion used  to  define  the  two  kinds  of  variation  treated:  random  and  sys- 
tematic variation.  If  the  environments  are  considered  to  be  defined  by 
first  order  (mean)  and  second  order  (variance)  moments  of  distributions 
which  are  better  represented  in  log -space  than  in  linear  space,  then 
the  log-normal  distribution  is  both  physically  plausible  and  convenient. 
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Environments  may  include  both  ambient  and  in  situ  conditions,  as 
well  as  characteristics  of  the  hostile  environment.  For  example, 
geologic  properties  such  as  depth  to  water  table  may  be  treated  as  an 
environment.  The  water  table  depth  could  be  a very  significant  param- 
eter for  components  sensitive  to  facility  flooding.  Hostile  environ- 
ments such  as  peak  values  of  free  field  displacements,  velocities  and 
accelerations,  or  shock  spectra  parameters  when  defined  at  appro- 
priate depths,  can  be  related  to  critical  internal  shock  environments 
or  system  responses.  Ground  shock  loads  on  protective  structures 
may  be  related  to  free  field  stresses  and  so  on.  In  general,  environ- 
ments are  a mix  of  ambient  and  in  situ  conditions  transformed  into 
hostile  environments  by  the  energy  associated  with  a catastrophic 
event. 

Among  the  most  important  aspects  of  environment  definition  is 
the  derivation  of  covariance  matrices  for  random  and  systematic 
variations  on  environment  predictions.  Covariance  or  correlation  is 
important  because  non-zero  covariance  between  two  environments  can 
strongly  affect  the  joint  probability  that  both  environments  exceed 
some  set  values  at  the  same  time.  For  a series  system  where  two 
components  are  affected  by  strongly  correlated  environments  (e.  g.  , 
debris  and  crater-related  ground  shock),  then  either  both  components 
fail  or  both  survive,  and  the  probability  of  system  kill  does  not  com- 
pound itself  as  more  components  are  considered.  Alternatively,  if 
a series  system  is  subjected  to  anticorrelated  hostile  environments 
(e.  g.  , air  induced  velocity  and  crater-related  velocity),  then  the 
probability  of  system  kill  can  be  substantially  greater  than  if  the 
environments  were  uncorrelated.  The  effects  of  covariance  on  parallel 
and  series-parallel  systems  are  also  important. 

Noticing  the  importance  of  covariance,  how  does  one  go  about 
deriving  covariance  matrices  for  random  and  systematic  variations  on 
the  environment  predictions?  The  answer  depends  upon  the  problem 
solving  mode. 

For  a well  known  system  in-place  at  many  sites,  one  can  perform 
environment  parameter  predictions  based  upon  known  in  situ  properties 
and  computer  code  calculation  results  and/or  empirical  scaling  laws. 
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A vector  of  pertinent  environments  can  be  constructed  for  each  site, 
and  multivariate  log-normal  statistics  can  be  computed  (Equation  (2-6)) 
for  one  or  more  groups  of  these  sites  to  define  the  mean  vector(s)  and 
random  covariance(s)  of  the  environment  predictions.  Next,  one  exam- 
ines the  environment  prediction  equations  (which  are  functions  of  in 
situ  conditions  and  the  manner  of  energy  application)  and  determines 
the  underlying  systematic  variation  or  uncertainty  on  material  proper- 
ties, loads  and  analysis  (prediction)  technique.  If  the  prediction 
equations  are  linear  (in  log -space),  the  environment  systematic 
variation  covariance,  £uu,  may  be  computed  directly  from  the 
covariance  of  underlying  variables,  £ 


vv 


uu 


L £ L 
vv 


(2-19) 


where  L is  the  matrix  which  defines  the  linear  prediction  technique  in 
log -space. 


For  example,  consider  a system  sensitive  to  only  two  environ- 
ments; vertical  velocity  (vz)  and  horizontal  velocity  (vh)  defined  by 


P P 

P c’  vh  ” P U 
m m 


where  P is  the  airblast  pressure,  U is  the  air  shock  speed,  p is  the 
media  density  and  c is  the  media  stress  wave  speed.  In  log-space  the 
prediction  equations  are 


log  vz  = log  P - log  pm  - log  c 


log  vh  = log  P - log  pm  - log  U 


Thus, 
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where  the  subscripts  1,2,  3,4,  refer  to  log  P,  log  p , log  c and  log  U, 
respectively. 

An  alternate  mode  of  generating  random  covariance  matrices  for 
environments  is  to  estimate  the  random  variation  of  underlying  variables, 
and  propagate  these  random  covariances  through  the  prediction  equa- 
tions as  was  described  for  systematic  variations. 


Scaling  laws  for  nuclear  weapons  effects  have  been  included  in 
FAST,  in  the  forms  specified  in  Equations  (2-7a  and  2-7b).  For  sur- 
face bursts,  these  forms  of  scaling  are  adequate  for  reasonable  inter- 
vals of  range,  R (and  corresponding  overpressure,  P);  yield,  W;  and 
height  of  burst,  H.  The  forms  of  scaling  have  the  property  that  logarith- 
mic covariances  are  constant  while  logarithmic  means  vary  linearly  as 
the  logs  of  overpressure,  range,  or  yield.  More  complex  forms  of  scal- 
ing both  the  mean  and  the  covariance  are  possible.  However,  only  simple 
scaling  laws  have  been  implemented  in  the  FAST  code  because  they  have 
proved  to  be  adequate  in  all  past  applications.  The  effects  of  complicated 
scaling  can  be  modeled  by  computing  as  many  baseline  threat  mean  vec- 
tors and  covariance  matrices  as  required,  or,  should  some  new  applica- 
tion dictate,  the  scaling  can  be  changed  as  it  has  been  for  nonlinear 
overpressure  scaling. 

Scaling  for  height  of  burst  is  accomplished  by  scaling  pressure 
and  related  effects  with  overpressure  and  yield  separately  for  each  H, 
and  by  scaling  crater-related  effects  as  a function  of  H,  based  on  the 
crater  volume  scaling  V = V (H). 
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2.  3.  2 Transfer  Functions 

Transfer  functions  normally  take  a single  free  field  environment 
and  convert  it  into  a critical  local  response.  To  illustrate  a more 
complex  case,  consider  a buried  structure  which  is  subjected  to  verti- 
cal overpressure,  P,  on  the  roof  and  horizontal  pressure,  KQ  P,  on 
the  walls.  The  lateral  pressure  coefficient  KQ  is  a property  of  the  in 
situ  soils,  and  as  such  proper  correlation  of  environments  requires 
that  K0  he  treated  as  an  environment  although  it  obviously  would  not 
scale  with  pressure,  yield,  or  HOB.  Now  the  fragility  can  be  expressed 
as  a function  of  the  dimensionless  static  ratio  of  stress  to  allowable 
stress,  S, 


Sr  = T- 

r 


for  the  roof 


(2-20a) 


s = tt- 
w V 


* K P 
w o 


for  the  wall 


(2  -20b) 


where  X and  X.  denote  the  static  stress  concentration  factors  for  the 
r w 

roof  and  the  wall  and  V and  V denote  the  allowable  shear  strengths 

r w 

of  the  concrete  roof  slab  and  the  concrete  wall.  Based  on  construction 
data,  one  might  find  or  anticipate  that  and  are  correlated  due 
to  workmanship,  weather,  or  other  factors.  In  order  to  preserve 

this  correlation  one  could  define  P/V  and  PK  /V  as  environments. 

' r o'  w 

Carrying  the  correlations  even  further,  the  stress  concentration 

factors  X.  and  X.  might  have  been  derived  from  the  same  set  of  data 
r w 

on  shear  failures  of  deep  slabs.  In  this  case,  an  error  in  data  inter- 
pretation would  manifest  itself  in  both  X.'s,  hence  these  correlations 
might  be  best  preserved  by  calling  X.rP/Vr  and  X^I^P/V^  environments. 
Then  the  transfer  function  defines  only  the  dynamic  amplification  fac- 
tors for  the  roof  and  wall,  respectively,  and  the  fragility  curves  might 
be  very  simple  curves  centered  on  the  dimensionless  constant  1. 
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It  is  only  the  ratios  of  response  to  allowables  which  matter. 
Statistics  on  these  ratios  will  be  reasonable  if  the  steps  which  are  used 
in  developing  environments,  transfer  functions,  and  fragilities  are 
reasonable.  Transfer  function  effects  should  not  be  arbitrarily  lumped 
into  either  environments  or  fragilities.  Nevertheless,  the  distinctions 
between  environments,  transfer  functions,  and  fragilities  are  sometimes 
left  to  the  convenience  of  the  user,  and  one  should  not  let  semantics 
limit  one's  flexibility  in  deriving  inputs  for  the  FAST  code. 


2.  3.  3 Fragilities 

Derivation  of  fragilities  starts  with  the  frequently  difficult  task  of 
identifying  the  critical  response  for  each  component.  Fragility  curves 
are  then  determined,  using  analysis,  test  data  or  a combination  of 
analysis  and  test.  These  fragility  curves  represent  the  best  estimate 
of  the  actual  failure  probability.  Each  fragility  must  be  carefully 
developed  to  ensure  that  it  is  not  too  low  or  too  high.  In  general, 
engineers  and  analysts  tend  to  produce  fragilities  which  are  conservative 
(too  low).  Special  precautions  are  strongly  recommended  to  provide 
fragilities  which  are  realistic.  Specific  examples  provide  the  best  guide 
to  understanding  fragility  derivation,  and  two  examples  are  highlighted 
in  the  sample  problem  of  Section  6.4. 

In  FAST  there  is  considerable  flexibility  in  the  representation  of 
fragilities.  This  flexibility  facilitates  specialization  where  component 
or  subsystem  fragilities  are  derived  by  experts  independent  of  the 
derivation  of  other  fragilities,  transfer  functions,  and  environments. 
However,  trouble  can  ensue  if  a transfer  function  is  misidentified  as  a 
fragility.  When  the  problem  is  properly  subdivided  and  environments 
and  transfer  functions  are  properly  derived,  the  development  of 
fragilities  is  facilitated. 

Azimuth  sensitivities  can  be  introduced  in  series  chainwise 
fashion,  which  enables  use  of  multiple  azimuth  tables  to  model  complex 
situations.  One  source  of  significant  azimuthal  effects  are  drag  loads 
on  blast  doors  or  non-surface  flush  and  non-axisymmetric  structure 
conditions.  Non-flush  surfaces  may  contain  multiple  failure  modes 
(i.  e.  , blast,  thermal,  debris,  shock,  radiation)  and  all  these  modes 
usually  occur  concurrently. 
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2.3.4  Network 

Development  of  the  system  network  involves  determination  of  the 
series /parallel  relationship  among  component  fragilities.  In  addition, 
the  manner  in  which  components  are  grouped  to  form  subsystems  needs 
to  be  established.  A few  of  the  possible  approaches  to  organizing  a 
study  are  to  group  into  logical  subsystems 

a.  components  which  perform  related  functional  relationships, 

b.  components  which  are  affected  by  related  environments, 

c.  components  which  are  in  nearby  locations,  and 

d.  components  which  have  similar  attack-failure-repair 
c hronology. 

The  manner  in  which  components  are  grouped  can  be  quite 
important  since  the  survivability  statistics  on  subsystems,  as  well  as 
on  components  and  on  the  system,  are  output  from  FAST.  By  the  use 
of  ganging  it  is  possible  to  run  a problem  in  which  one  system  is 
organized  from  several  logical  viewpoints.  Often  this  combined  approach 
will  maximize  the  user's  understanding  of  the  survivability  statistics. 

In  order  to  model  component  fragility  over  a wide  range  of 
conditions,  it  is  often  necessary  to  expand  the  system  network.  For 
example,  a component  may  be  sensitive  to  displacements  for  high 
weapon  yield  conditions  and  velocities  for  low  weapon  yield  conditions. 
Such  a component  would  be  best  modeled  by  two  fragilities  in  a series 
or  parallel  relationship  (depending  on  the  specifics  of  the  example). 
Impacts  between  objects  and  combinations  of  stress  and  deformation 
also  can  result  in  complicated  logic.  Sometimes  network  logic  complica- 
tions can  be  removed  or  reduced  by  using  more  complex  transfer 
functions  and/or  more  complex  fragilities.  The  requirements  of  the 
study  and  characteristics  of  the  system  should  play  an  important  role 
in  deciding  how  the  modeling  is  performed. 

Other  changes  from  the  functional  system  network  to  the  FAST 
system  network  are  dictated  more  directly  by  the  questions  which  the 
FAST  calculations  are  expected  to  help  answer.  For  example,  one 
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might  want  to  find  out  how  often  Component  A and  Component  B both 
fail  even  though  the  failure  of  either  would  result  in  system  failure. 

To  calculate  joint  failure  probabilities,  one  could  construct  a dummy 
system  (using  the  ganging  technique)  containing  Components  A and  B 
in  parallel. 

In  finalizing  a network,  one  should  prune  a preliminary  fault 
tree,  making  use  of  quantitative  knowledge  of  fragilities  and  environ- 
ments to  eliminate  superhard  components  and  components  which  almost 
always  fail  in  conjunction  with  other  more  critical  components.  In  this 
manner  complex  systems  are  meaningfully  reduced  to  manageable 
proportions.  Complete  identification  and  documentation  of  even 
superhard  items  should  be  accomplished  as  part  of  the  network  synthesis 
prior  to  pruning  the  fault  tree. 
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3.  0 INPUT  DESCRIPTION 


This  section  defines  all  of  the  inputf.  required  for  the  FAST  Code 
and  describes  the  required  input  format. 

The  inputs  are  contained  on  five  classes  of  cards 

• Fragility  - Transfer  Function  Cards 

• System  Network  Cards 

• Environment  Parameter  Cards 

• Threat  Condition  and  Convergence  Criteria  Cards 

• Miscellaneous  and  Control  Cards 

Each  kind  of  input  card  will  be  discussed,  giving  the  card  format 
and  a table  defining  all  variables.  Each  format  is  illustrated  by  cards 
taken  from  the  Sample  Problem  baseline  run  in  Appendix  A,  unless 
otherwise  indicated. 

In  discussing  the  inputs,  any  number  that  is  not  identified  as  an 
integer  will  be  a real  (floating  point)  number.  All  integers  should  be 
written  without  a decimal  point  or  exponent.  Real  numbers  may  be 
written  in  any  of  the  equivalent  forms. 


i 

I 


i 


1 Z 3 123.  123.  0 = 1.  23E2 

or  -3.  2E-1  -0.32 

3.  1 FRAGILITY  - TRANSFER  FUNCTION  CARDS 

The  basic  parameters  that  can  be  defined  on  the  fragility  - transfer 
function  card  are  summarized  in  Table  3-A.  For  convenience,  this 
card  is  referred  to  as  a component  card.  The  basic  form  of  the  com- 
ponent card  is 


Ci  /0’  fPL  i’V  fPu  /p f2>  fPL>2’fp2’  fpU,2/ 

* ,1  * 

...  /I,  i > f , f /RK  - ne  , tf.,  tfu.  $ 

PL,  6 p6  PU,  6 1 1 1 


The  information  from  a typical  component  card  is  reproduced 
below  with  a description  of  each  parameter  on  the  card. 


C101  /0,  30,  35,  40/1,  tO,  45,  50/R1  = 13,  .2,  1.2$ 

where  C101  is  the  component  card  identification.  Components 
are  always  identified  by  integer  numbers. 
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Table  3 -A.  Fragility  - Transfer  Function  Parameters 


Parameter 


fPL,j’  V Vj 


Description 

Integer  designating  component  number 

The  probability  of  failure  at  the  level 

Piecewise  linear  representation  of  the 
fragility.  in  , f , f are  the  cor- 

j Pj  F 


are  the  cor- 


PL,j  Pj  PU,j 

responding  lower  bound,  nominal  and 
upper  bound  values  of  the  component 
response,  respectively,  at  pf..  A maxi- 
mum of  6 points  can  be  considered. 

f , and  f need  not  be  specified  if 
PL,j  PU,j 

there  is  no  systematic  variation 

Environment  parameter  reference  number 
corresponding  to  component  i 

Transfer  function 

Transfer  function  k factor 

Storage  register  for  transformed 
environment 

Terminates  component  card  and  may  be 
followed  by  comments 


*Neeo  not  be  specified  if  equal  to  unity 


/0,  30,  . 5,  40/1,  40,  45,  50/ represents  the  following  fragility  curve: 


C101 

l.Oi 


40  45  50 


/// 


30  35  40 

Component  Response 


ne  = 13 
tf  = 0. 2 
tfu  = 1.2 


i 


•13 


ft 


F 

I 


ii ....  i Ji  m mm,nimwmm.,mwfl  "*->Mff,,  »iV*n^r"Siwwr«w™"'ir-.-.^7' ■'•  ""mM'.'w.'u.* 'i'*.' .mm; >•■ 
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R1  = 13,  . 2,  1 .2  describes  the  relationship  between  an  environ- 
ment and  a component  response.  Environments,  like  components, 
are  referenced  by  integer  numbers.  In  this  car.e,  the  reference 


is  to  environment  number  13.  Of  Ej^  is  the  value  associated  with 


environment  13,  then  the  transfer  function,  tf  = 0.  2,  is  used  as 
follows  and  the  result  temporarily  stored  in  R1 


R 1 = E^  *0.  2 


The  transfer  function  k factor,  tfu  = 1.  2,  applies  bias  to  Rl.  The 
result  is  then  stored  back  into  Rl. 


The  information  on  a component  card  could  be  as  simple  as 
C604  = /0,  10/1,  10/R1  = 6,  1,  1 $ 


which  represents  a vertical  (cookie  cutter)  best  estimate  fragility.  The 
component  response  is  equal  to  the  unmodified  environment  6,  since  the 
tf  and  the  tfu  are  both  equal  to  unity.  It  is  important  to  note  that  all  fra- 
gilities must  be  defined  at  = 0.  0 and  = 1.  0.  The  fragility  curve  for 
this  component  is 


C604 


ne  = 6 

tf  = 1.0 

tfu  = 1 . 0 
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Rl 


Component  Response 


The  tf  and  tfu  values  need  not  be  punched  when  they  are  equal  to 
1,  so  the  following  forms  are  also  equivalent  for  C604: 


C604  = /0,  10/1,  10/R1  = 6 $ 

C604  = /0,  10/1,  10/Rl  = 6,  1 $ 

C604  = /0,  10/1,  10/Rl  = 6,  , 1 $ 


Fragility  curves  and  transfer  functions  that,  are  more  complex  can 
be  handled  in  FAST.  These  are  discussed  in  the  following  section. 

This  section  can  be  skipped  during  initial  acquaintance  with  the  code. 
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3, 1,  1 Special  Functional  Operations 


The  environment  field  can  have  more  complex  expressions 
including  the  use  of  CAN  and  USER.  CAN  provides  the  option  to  s.dd 
environments  or  the  transformed  environments  and  USER  provides 
azimuth  sensitivity.  The  first  example  illustrates  the  use  of  CAN 
(Appendix  A does  not  include  an  example): 


C502  = /0,  2,  5,  7/1,  3,  6,  ?/P-l  = 3/R2  = 6/R3  = 2,  0.  2/CAN(3)/ 
R1  = Rl,  1,  1.  2 $ 


where  Rl  = 3 stores  environment  3 in  Rl 


R2  = 6 stores  environment  6 in  R2 


R3  = 2.  , ,2  stores  0.  2 times  environment  2 in  R3 

CAN(3)  indicates  how  many  environments  to  add  and  stores  the 

result  in  Rl,  Thus,  Rl  = Rl  + R2  + R3. 

Rl  = Rl,  1,  1.2  applies  the  transfer  function  k factor  to  Rl  and 
stores  result  back  into  Rl. 


The  input  USER(IG)P)  supplies  the  calling  sequence  parameter  I0P 
to  a subroutine  called  USER  which  must  be  provided  by  the  program  user 
in  order  to  account  for  azimuth  sensitivity.  The  parameter  I0P  then 
allows  for  azimuth  table  selection.  An  example  of  such  a subroutine  is 
given  in  Appendix  A. 


Component  card  Cl 02  illustrates  a reference  to  subroutine  USER. 


C102  = /0,  9,  10,  11/1,  11,  12,  13/R1  = 17/USER  (1)=  Rl,  1.0,  1.  2$ 


where  Rl  = 17  indicates  environment  17  is  stored  in  Rl 


USER( I ) indicates  the  azimuth  table  to  apply  for  ne^  in  C102 
(see  the  following  for  example). 


Azimuth  table  taken  from  subroutine  USER  (see  Appendix  A): 


XI  /0,  30,  45,  56,  60,  70,  90,  110,  120,  124,  135,  150,  180,  210, 
225,  236,  240,  250,  270,  290,  300,  304,  315,  330,  360/ 


Y1  /l,  .866,  .707,  .552.  .578,  .627.  .667,  .627,  . 578,  .552,  .707,  .866, 

1,  .866,  .707,  .552,  . 578,  ,627,  .667,  .627,  . 578,  .552,  .707,  .866,  1.  / 


The  previous  table  represents  the  following: 


Azimuth 
Coefficient,  ■ 
AZ  (4>) 

(i.  e.  , Yl) 

0 L 


To  be  consistent  with  CAN,  USER  stores  the  result  in  R1 

Rl  = Rl,  1,  0,  1. 2 indicates  the  k factor  is  applied  and  the  result 

stored  in  Rl. 


3.  2 SYSTEM  NETWORK  CARDS 


Table  3-B  summarizes  the  Boolean  network  parameters. 


Table  3-B.  Network  Parameters 


Parameters 

Example 

Description 

+ 

Boolean  operator  for  combining 
series  elements 

* 

-Cg> 

Boolean  operator  for  combining 
parallel  elements 

( ) 

►0-0 

Used  to  order  how  elements  are 
combined 

Ci  or  i 

C101  or  101 

An  integer  that  identifies  a 
component  element.  The  C need  not 
be  punched 

SSj 

SS  1 

An  integer  than  identifies  a sub- 
system element 

SYSk 

SYS  100 

An  integer  that  identifies  a system 
element.  The  results,  when 
printed,  are  identified  by  negative 
k (-k) 

$ 

Terminates  the  equation 
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30  60 

Azimuth,  4> 


90  120  150 

(degrees)  (i.  e.  , XI) 


180 


E 


t % 


ipjpppppi  .u  !i. 


The  first  two  examples  are  selected  from  Appendix  A. 
SSI  = 101  + 102  * 103  $ 


The  operator  * takes  precedence  over  + , thus  the  probabilities  of 
failure  P^^  * ^io3  are  com^ne<^  first  to  compute  a combined  proba- 
bility of  failure,  Py.  Then  Pr  is  combined  with  Pjqj  as  follows: 


SSI 


= 1 - U-Pioi*  (1-p-> 


Parentheses  are  used  to  change  the  order  of  operations: 
SS2  = (131  + 132)  * 133  $ 


SD-DS 

1 133 


In  this  case  Pj^j  and  are  combined  first  to  compute  a com- 

bined probability  of  failure:  Pf  = 1 - (1_Pi3i)  ^"P132^  ThenPris  combined 


with  P233  as  follows: 


p _ p * p 

SS2  r 133 


In  each  of  these  examples  the  result  is  designated  as  a subsystem 
(SSI  and  SS2),  but  they  could  have  been  systems  (SYS1  and  SYS2).  The 
following  examples  are  not  typical,  but  they  illustrate  that  parentheses 
can  be  nested  (to  almost  any  depth)  and  any  elements  can  be  combined 
as  shown  by  the  following  equations: 


SS6  = SYS7  + SS10*((10+11)*8  + 16)  $ 
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SS6  = SYS7  + SS10*((10+  11)#8  +16)  $ 


SYS  10  = SYS8  + SS6  + 21  $ 


|SYS8| 


There  are  two  network  restrictions.  First,  just  as  components 
must  be  entered  before  they  are  used  in  an  equation,  subsystems  and 
systems  must  be  defined  (equations  must  be  in  logical  order  as  in  a 
FORTRAN  program).  Second,  redundant  parentheses  are  not 
permitted: 

(101)  $ 

((101  + 102))*103  $ 

(101  + 102  + 103)  $ 

3.  3 ENVIRONMENT  PARAMETER  CARDS 

Environments  are  entered  on  various  cards.  Each  environment 
parameter  is  described  in  Table  3-C. 

Each  of  the  vectors  in  Table  3-C  can  define  up  to  35  environments. 
The  only  card  that  is  not  illustrated  in  Appendix  A is  the  ENV  card. 

This  card  is  an  alternate  method  of  entering  and/or  changing  environ- 
ment parameters;  the  general  format  is: 

ENV  i r MODE(i),  ULN(i),  K-FAC(i),  ALPHA(i),  BETA(i), 
GAMMA(i),  HB(i),  GAMMAB(i) 

(The  last  three  parameters  need  not  be  entered  if  they  are  not  used; 
just  as  GAMMAB  and  HB-GAMMAB  are  not  required.  ) 

To  change  originally  punched  and  entered  environment  param- 
eters, the  following  cards,  which  supersede  the  previous  information, 
are  entered  after  the  original  cards 

ENV2  = 1,  . 5,  2,  . 1,  . 05,  .8,  . 02,  .1 
ENV 3 = 2.  .6.  1.5.  .2.  .1.  .7.  .2.  .7 
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Table  3-C,  Environment  Parameters 


Parameter 


MODE 


Description 

Vector  which  identifies,  by  integer  numbers,  the 
type  of  environment  and  scaling  option 


MODE  Scaling  Option  Type  of  Environment 


ULN 


GAMMA 


range 

range 

range 

overpressure 

overpressure 


log-normal 

normal 

normal  (log  base  10) 

log-normal 

normal 


6 overpressure  normal  (log  base  10) 

Vector  of  logarithmic  environment  medians 

Vector  of  yield  scaling  exponents 

Vector  of  range  or  pressure  scaling  exponents 

Vector  of  coefficients  to  scale  height  of  burst 
(HOB/Wr/3) 


HB-GAMMAB  Height  of  burst  break  point  and  scaling  coefficient 
are  entered  in  pairs  into  a vector 

S Used  to  enter  rows  or  columns  of  symmetric 

random  variation  covariance  matrix 


K-FAC 


Vector  of  K factors  for  systematic  environment 
variation 


Used  to  enter  rows  or  columns  of  symmetric 
correlation  matrix 


ENV 


Alternate  parameter  for  entering  several 
environment  parameters 


or  the  equivalent  alteration  cards  are  entered  after  the  original  cards 


MODE(2)  = 1,  2,  ULN(2)  = .5,  .6,  K-FAC(2)  = 2,  1.5 
ALPHA<2)  = . 1,  . 2,  £ETA(2)  = . 05,  . 1,  GAMMA(2)  = . 8,  . 7 
HB -GAMMA B(2)  = . 02,  .1,  .2,  .7 

As  Appendix  A shows,  the  formats  for  the  environment  param- 
eters are  not  complex.  Thus,  the  MODE  card  will  be  discussed  to 
represent  the  vector  inputs  and  the  S cai'd  will  be  discussed  to  represent 
the  matrix  inputs. 


MODE  = 16  * 4,  3*1 

This  is  a convenient  format  that  is  equivalent  to: 

MODE  = 4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  1,  1,  1, 

The  integer  number  following  the  S identifies  the  diagonal  element 
of  the  matrix.  For  example: 

Si7  - . 0230,  0.  107,  . 0180 
which  defines  the  following  elements 

s17  17  = 
s18, 17  = °*  107 
s 1 9, 17  = ,018° 

3.4  THREAT  CONDITIONS  AND  CONVERGENCE  CRITERIA  CARDS 

Table  3-D  defines  the  threat  condition  parameters. 

The  choice  for  the  sample  problem  was  to  enter  the  following 
cards  to  describe  the  threat. 

NOMFWH  = 6 00,  5 $ 

where  P = 600,  W =5  and  H = 0 
o o o 


and 


W = 5 


OP  = <;00,  400,  500,  600,  700,  800,  1000,  1200 
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Table  3-D,  Threat  Condition  Parameters 


Parameter 


Description 


NOMRWH 


Name  used  to  identify  the  input  parameters 

R , W and  H . NOMRWH  = R , W , H 
0 0 o o’  o*  o 


NOMPWH 


Name  used  to  identify  the  input  parameters 
P0,  WQ  and  H0.  NOMPWH  = P0  W , H . 

Only  one  of  these  two  cards  (NOMRWH  or 
NOMPWH)  is  required 


Nominal  range 


Nominal  yield 


Nominal  height  of  burst  - assumed  to  be  zero 
if  not  entered 


Nominal  overpressure 


Ranges  of  interest  (see  OP) 


Yield  of  interest 


Height  of  burst  of  interest  - assumed  to  be 
zero  if  not  entered 


Vector  of  overpressures  of  interest.  Either 
R or  OP  is  given,  not  both 


for  yield  of  5 MT  and  various  pressures  (ranges  are  calculated  auto- 


matically by  the  program  as  is  R^). 


Table  3-E  defines  the  convergence  criteria  parameters. 


For  the  sample  problem,  the  choice  was  to  use  the  default  values. 
The  default  values  are  the  values  automatically  used  by  FAST  if  no 
other  specific  values  are  included  in  the  input.  The  default  values,  in 
the  same  order  presented  in  Table  3-E,  are  the  same  as  entering  the 
following  card. 


CONFI  = 5,  50,  50,  . 05,  . 03,  . 02,  .5 


_•  4 ..’all ?— 


■ 


Table  3-E.  Convergence  Criteria  Parameters 


3.  5 CONTROL  CARDS  AND  MISCELLANEOUS  CARDS 

The  control  cards  cause  the  program  to  perform  certain  opera- 
tions, which  are  completely  described  in  Table  3-F, 

Table  3-F.  Control  Cards 


tw  .mmrnm 


Card 

Description 

RUN$ 

Indicates  that  all  the  information  for  a case  has  been 
entered. 

RESETS 

Indicates  the  beginning  of  all  new  information  — pre- 
vious information  is  deleted. 

STOPS 

Indicates  the  end  of  data  and  causes  the  program  to 
stop. 

PLOTON-  s$ 

Turns  the  printer  plot  option  on,  where  s is  the 
overpressure  per  inch  to  scale  the  plot  abscissa, 
which  is  six  inches  long. 

PLOTS 

Indicates  when  to  plot  the  results  Lorn  one  or 
more  cases. 

PLOTOFFS 

Turns  the  plot  option  off. 

Parameter 

Description 

CONFI 

Name  of  the  card  which  provides  the  parameters 
CONFI=  NIMIN,  NIMAX,  NOMIN,  EPSI,  EPSO, 
EPSOC,  FRACTILE 

NIMIN 

Number  of  inner  loop  iterations  to  do  before  test- 
ing inner  loop  convergence  (integer  number) 

NIMAX 

Number  of  outer  loop  it  -rations  to  do  before 
accepting  the  inner  loop  convergence  (integer 
number) 

NOMIN 

Number  of  outer  loop  iterations  before  testing 
outer  loop  convergence  (integer  number) 

EPSI 

Required  standard  deviation  for  inner  loop 
convergence  test 

EPSO 

Required  standard  deviation  for  convergence  of 
outer  loop  mean  (first  outer  loop  test) 

EPSOC 

Required  standard  deviation  for  convergence  of 
median  (0.  5)  or  more  fractile  (second  outer  loop 
test) 

FRACTILE 

Vector  of  fractiles  from  1 through  7 

1 P'M  M.T  fi  "1— 
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Table  3-G  defines  some  miscellaneous  cards  and  features.  Some 
of  these  cards  are  not  commonly  used  but  are  available. 


Table  3-G.  Miscellaneous  Cards  and  Features 


Card/ 

Feature 

Description 

RANUN 

RANUN=NUMUR,  NUMNR,  On  some  computer 
systems  the  starting  values  for  the  random  num- 
ber generators  (uniform,  and  normal)  can  be 
changed,  but  the  user  must  be  familiar  with  FAST 
and  the  available  generators.  The  default  values 
assumed  by  the  TRW  random  number  generators 
are  1. 

RECODE 

This  card  allows  the  data  to  be  punched  on  both 
an  029  and  026  keypunch  machine 

Ci=$ 

Deletes  Component  i 

SSj=$ 

Deletes  Subsystem  equation  j 

SYSh=$ 

Deletes  System  equation  h 

null -field 

A null-field  is  described  below 

L format 

For  MODE=3  or  6 the  ULN  fragility  anu  transfer 
function  will  be  in  log  base  10.  If  the  L format  is 
used  in  place  of  the  E,  the  program  will  take 
lo8l0- 

18L0  = 1.8  LI  s 180L-1 
where  the  result  is: 

1.  8Ll  = log10  (1.  8*1 01 ) = 1.  25527 

The  RECODE  card  is  fixed  format  (no  blanks)  that  starts  in 
column  1. 


RECODE=ABCDEFGHIJKLMNOPQRSTUVWXYZ01 23456  789,+  -x/.  ()^$ 


If  the  data  deck  has  been  punched  on  both  an  029  and  026  keypunch, 
then  the  RECODE  card  will  be  useful.  In  front  of  029  cards,  place  a 
RECODE  card  punched  in  029  and  an  026  RECODE  card  in  front  of 
026  cards. 
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A null-field  is  defined  as  a numeric  field  that  does  not  contain  a 
number.  This  is  useful  for  changing  values  on  cards  like  CONFI  or 
ENV.  As  an  example,  this  card  will  only  change  the  values  for  NOMIN 
(3rd)  and  EPSO  (5th) 


CONFI=2*,  40,,  . 02 


3.  6 CASE  STACKING 

The  FAST  inputs  can  be  manipulated  to  do  ganging,  job  stacking 
and  case  stacking.  Ganging  permits  several  system  configurations  in 
one  case.  The  RESLT$  card  deletes  previous  information,  which  per- 
mits stacking  jobs  in  a computer  run. 

Appendix  A (yield  and  HOB  perturbations)  shows  a typical  example 
of  case  stacking,  where  the  yield  and  height  of  burst  are  varied.  In  other 
words  the  information  in  a previous  case  is  available  for  the  next  case  and 
can  be  modified  as  required.  For  example,  components  and  network 
equations  can  be  replaced  or  deleted.  Thus  C604,  from  Appendix  A, 
(Baseline  System  Sample  Problem)  could  be  changed  to 

C604=/0,  10/1,  10/R1  = 6,  1.2$ 
or  deleted  (be  sure  to  modify  network)  as  shown  below: 


C604=$ 
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4.0  OUTPUT  DESCRIPTION 

All  of  the  input  data  is  printed  back  out.  Since  the  "USER(IOP)" 
subroutine  is  an  integral  part  of  the  input,  the  subroutine  is  printed  out 
to  show  the  azimuth  data;  input  cards  are  printed  to  show  information  on 
the  component,  network,  and  environment  cards. 

If  there  are  any  keypunch  errors  or  errors  of  omission  or  incon- 
sistencies, error  messages  are  printed  out.  See  Section  4.  2 for  a list 
of  the  messages  and  suggestions  for  correction.  Finally,  the  actual 
calculated  data  for  the  probability  of  survival  for  each  component,  sub- 
system, and  systems  are  displayed. 

4.  1 SURVIVABILITY  RESULTS 

If  the  rank  for  either  the  random  (S)  or  systematic  (A)  variation 
matrix  is  less  than  the  dimension,  then  the  following  illustrates  what  will 
be  printed. 

S MATRIX  - DIMENSION  12  HAS  RANK  11 

For  each  range  and  overpressure,  the  input  ULN  vector  is  scaled 
and  the  median  values  of  the  environment  are  displayed  as  follows: 


ULN(  0) 
ULNI  1) 
ULN I 2 ) 
ULNI  4) 
ULNI  51 
ULNI  61 
ULNI  8) 
ULNI 10) 
ULNI 11) 
TJLIU27 
ULNI 13) 

ULNI 15) 
ULNI  6) 
ULNI 17 » 
ULNI It  I 
ULNI 19) 


8*  00 000E+02 
1.96435E+01 
5«  06580E+02 
7« A0276E+02 
1.53933E+01 
5.09239E»00 
7.73887E+01 
4.02720E»03 
2.07585E+03 
1.333UE+02 
2.64582E+02 
1.Q0056E+02 
8.  00056E+02 
8.00056E+02 
2. 18480E+01 
3. 56735E+01 
1. 79089E+01 


Note  that  ULN  (0)  is  not  an  input  uut  is  the  current  value  for  Overpressure, 
which  can  be  referenced  on  the  component  cards  as  environment  0,  i.e. 

R 1 0.  The  information  on  the  condition  of  the  threat  at  which  the  above 

environments  are  calculated  are  displayed  as  follows: 


where 


W-  5.  0000E+  00  H-  0.  R-  5.  1V09E-01  P=  6.  0000E+02 


W - weapon  yield,  MT 
H = height  of  burst,  K-ft 
R - miss  distance,  NM 
P - overpressure,  psi 


During  the  Monte  Carlo  process,  the  code  tests  for  convergence  as 

described  in  Section  2.  2.  5.  1,  If  there  is  no  convergence,  a line  of 

information  is  printed  which  includes  the  number  of  inner  loop  iterations 

(NI),  number  of  outer  loop  iterations  (NO)  and  the  test  criteria  (EPS). 

2 2 

EPS  is  (equation  2-1  3),  fQl  (equ  tion  2 - M ) or  2 € (equation  2- 1 5). 

If  the  process  converges  within  025  outer  loop  iterations,  then  a line 
is  printed  to  display  the  number  of  iterations.  A sample  appears  below. 
This  process  can  be  controlled  by  the  CONFI  input  card  (see  Section  3.  4). 


W* 

5*  0000E+00  H* 

0. 

R= 

4.6795E-01 

P*  8.0000E  *02 

NO 

CONVERGENCE  - 

NI* 

5 N0= 

0 EPS* 

5.0E-02 

NO 

CONVERGENCE  - 

NI* 

5 NO* 

i EPS* 

5.0E-02 

NO 

CONVERGENCE  - 

NI* 

5 NO* 

2 EPS* 

5.0E-02 

NO 

CONVERGENCE  - 

N I = 

5 NO* 

3 EPS* 

5.0E-02 

NO  CONVERGENCE 
NO  CONVERGENCE 
NO  CONVERGENCE 


^ HQr.  . 4 ?P5=  -5.QE-QZ. 

5 NU=  5 E PS=  5.0E-02 

5 NU=  6 EPS=  5.0E-02 

5 N0=  7 EPS=  5.0E-02 


NO  CONVERGENCE 
NO  CONVERGENCE 
NO  CONVERGENCE 
NO  CONVERGENCE 
NO  CONVERGENCE 

-NO  CONVERGER C£ 
NO  CONVERGENCE 
NO  CONVERGENCE 
NO  CONVERGENCE 
CONVERGENCE 


5 NO3 
5 N0= 
5 NQ  = 
5 NO3 
5 N0  = 

-5  NO 
5 NO 
5 NO 
5 NO 
5 NO 


41  EPS3 

42  EPS3 

43  EPS3 

44  EPS3 

45  EPS3 

46_E1>£. 

47  EPS 

48  EPS 

49  EPS 

50 


5.0E-02 

5.0E-02 

J5.0E-02 

5.0E-02 

5.0E-02 

5^QE--Q2 

5.0E-02 

5.0E-02 

5.0E-02 
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The  next  FAST  display  is  the  histogram  shown  in  Table  4-A.  The 
column  to  the  far  right  lists  the  input  numbers  which  identify  components  ^ 
systems  (negative  numbers)  and  subsystems.  There  are  twenty-two 
columns  to  store  probability  numbers  for  the  systematic  variation  distri- 
bution. The  first  column  is  for  the  count  of  the  probability  of  survival 
(P  ) equal  to  0 and  the  last  column  is  for  P =1.0,  while  the  twenty 
divisions  between  them  are  incremented  by  steps  of  0.  05  probability. 

The  number  of  outside  loops  determine  the  total  number  of  marks  that 
are  put  into  all  increments.  The  average  probability  of  the  inside  loop 
determines  which  increment  gets  the  mark  for  each  outside  loop  iteration. 

The  following  printout  section  gives  the  cumulative  distribution  of 
systematic  variation  by  individual  components,  subsystems  and  systems, 
and  are  ranked  by  the  50  percentile  or  median  of  the  distribution.  The 
component,  subsystem  and  system  names  (numbers)  are  listed  under 
NAME  on  the  left  hand  side  with  columns  d-j  ipnated  1,  10,  30,  50,  70, 

90,  99  for  the  percentile  of  the  distribution  for  each  component,  etc. 

An  interpolation  routine  is  used  which  prints  this  information  by  normal- 
izing to  the  number  of  ite-  .tions  actually  performed. 

individual  components 

NAME i 10 30 50 70 90 99_ 

133  0.000  0.000  0.000  0.000  0.000  0.000  1.000 

101  0.000  0.000  0.000  0.000  0.000  .200  1.000 


I Mi 

9 

11  MUM 

.82  5 

1.000 

1 

1.000 

— A a 

1.000 

AJ  iJMM 

1.000 

1.000 

k LUUU 

1.000 

j 

705 

0.  000 

0.000 

.850 

1.000 

1.000 

1.000 

1.000 

; 

INDIVIDUAL 

SUBSVST 

em£ 

NAME 

1 

10 

30 

50 

70 

90 

99 

1 

0.000 

0.000 

0.000 

0.000 

.400 

1.000 

1.000 

8 

0.000 

0.000 

0.000 

. 100 

.475 

.962 

1.000 

7 

.002 

.017 

.050 

. 167 

.340 

.583 

.92  5 

2 

0.  000 

0.000 

.250 

.650 

1.000 

1.000 

1.000 

5 

.275 

• 420 

.700 

.850 

1.000 

1.000 

1.000 

4 

0.000 

0.000 

• 850 

1.000 

1.000 

1.000 

1.000 

• 

6 

.825 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

3 

0.000 

0.000 

.850 

1.000 

1.000 

1.000 

1.000 

INDIVIDUAL 

SYSTEMS 

NAME 

1 

10 

30 

50 

70 

90 

99 

j 

-100 

0.000 

0.000 

0.000 

• 011 

.029 

.046 

.375 

1 

-1 

0.000 

0.000 

0.000 

.015 

.035 

.117 

.375 
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The  process  described  on  page  4-3  is  repeated  for  each  overpressure 
(OP)  or  miss  distance  (R)  designated  in  input. 

For  each  system  and  subsystem  designated  (components  are  not 
plotted),  a printer  plot  is  displayed  if  "PLOTON"  and  "PLOT$"  were 
input.  Again,  the  plots  are  identified  by  the  "minus"  number  for  system 
and  non-signed  number  for  subsystem.  The  abscissa  is  the  overpressure. 
The  ordinate  is  the  probability  of  survival  from  0 to  1.0.  The  points  are 
designated  at  the  vertex  of  the  letters  "A"  and  "V".  The  following  example 
shows  subsystem  8 plotted. 
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4.  2 ERROR  MESSAGES 

The  program  does  extensive  input  checking  for  errors  or  inconsis- 
tencies before  allowing  the  job  to  go  into  execution.  The  various  error 
mes;  ages  produced  along  with  an  identification  of  the  subroutine  from 
which  they  are  printed,  and  some  reasons  for  their  occurrence  are  listed 
below. 


PROGRAM  COULD  NOT  IDENTIFY  VARIABLE  NAME. 
FATAL  ERROR  PRECEDES  COLUMN  XXX  CARD  XXX. 
ERROR  SCAN  CONTINUES  WITH  EXECUTION  SUPPRESSED. 


Printed  from: 


NLNAME 


Explanation:  Look  for  (1)  missing  commas  or  (2)  other  format  errors,  or 
(3)  misspelled  word  (name). 


FATAL  ERROR  - EXPECTED  INTEGER  NUMBER  PRECEDING 
COLUMN  XXX  CARD  XXX 


Printed  from: 


NLINTG 


Explanation:  Look  for  (1)  missing  commas  or  (2)  non-integer  number. 


COULD  NOT  RECOGNIZE  LOCAL  ENV.  RELATION. 

FATAL  ERROR  PRECEDES  COLUMN  XXX  CARD  XXX. 
ERROR  SCAN  CONTINUES  WITH  EXECUTION  SUPPRESSED. 


Printed  from: 


COMPNT 


Explanation:  The  instructions  for  finding  component  response  from 
environment  came  across  a format  error  - usually,  a 
mispunched  "Ri",  "CAN"  or  "USER". 

TOO  MANY  ENV.  RELATIONSHIPS  FOR  INSC  TABLE. 
FATAL  ERROR  PRECEDES  COLUMN  XXX  CARD  XXX. 
ERROR  SCAN  CONTINUES  WITH  EXECUTION  SUPPRESSED. 


I 


% 


l 


wr"v.‘  * mv,  n.t.  miy  j~, 


Printed  from: 


COMPNT 


Explanation:  Too  many  calls  for  component  response  operations  - 
redimension  INSC. 


EQUATION  TOO  LONG  FOR  KEY  - DIVIDE  EQUATION. 
FATAL  ERROR  PRECEDES  COLUMN  XXX  CARD  XXX, 
ERROR  SCAN  CONTINUES  WITH  EXECUTION  SUPPRESSED 


Printed  from: 


BOOLE 


Explanation:  Network  expression  too  long  - watch  for  extraneous 
parentheses.  Note:  ’^operation  haj  precedence  over 
+ operation. 


NAME  NOT  FOUND  OR  SYMBOL  ERROR. 

FATAL  ERROR  PRECEDES  COLUMN  XXX  CARD  XXX. 
ERROR  SCAN  CONTINUES  WITH  EXECUTION  SUPPRESSED. 


Printed  from: 


BOOLE 


Explanation:  Watch  for  (1)  mispunch  or  (2)  missing  component  card, 
subsystem,  or  system  network. 


DIMENSION  PS  ARRAY  LARGER  FOR  BOOLEAN  SOLUTION. 
FATAL  ERROR  PRECEDES  COLUMN  XXX  CARD  XXX. 
ERROR  SCAN  CONTINUES  WITH  EXECUTION  SUPPRESSED 


Printed  from: 


LOGIC 


Explanation:  Network  has  too  many  parentheses  - redimension  PS 
array.  Note:  Do  not  place  parentheses  around  single 
term. 


TOO  MANY  INSTRUCTIONS  FOR  INSNET. 
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FATAL  ERROR  PRECEDES  COLUMN  XXX  CARD  XXX. 
ERROR  SCAN  CONTINUES  WITH  EXECUTION  SUPPRESSED. 


Printed  from:  LOGIC 

Explanation:  Too  many  BOOLEAN  OPERATIONS  (*,  + ) 


ENVIRONMENT  XXX  IS  MISSING.  EXECUTION  SUPPRESSED. 


Printed  from: 


SELECT 


Explanation:  At  least  the  "MODE"  for  the  environment  was  not  given 

the  user  should  check  all  environment  parameters. 


FATAL  ERROR  - 


jss  | 

(SYS) 


XXX  DELETED.  CHECK  NETWORK 


/SS  \ 
\ SYS) 


XXX 


Printed  from: 


NETCK 


Explanation:  Component,  subsystem  or  system  missing  but  referenced 

in  subsystem  or  system. 


WARNING  - COMPONENT  XXX  IS  NOT  REFERENCED. 


Printed  from: 


NETCK 


Explanation:  The  component  is  defined  by  a component  card  but  none  of 

the  network  equations  referenced  the  component. 

Note:  Non-fatal  error,  execution  proceeds. 


SYS  XXX  USES 


jss  j 

[SYS) 


XX  XXX  TIMES 


NETWORK 


PROBABILITY  IS  INCORRECT. 
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NETCK 


Printed  from 


Component,  subsystem,  or  system  is  used  more  than  once 
in  basic  system  expression. 


FATAL  ERROR  - HOB  TOO  LARGE 


Printed  from 


Height  of  burst  is  too  large  for  the  range  and/or  yield 
of  construction. 


PARENTHESES  DO  NOT  BALANCE 


Printed  from 


Explanation:  Self-explanatory 
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5.  0 COMPUTER  PROGRAM 

This  section  and  Appendix  B describe  the  FAST  computer 
program.  Appendix  B contains  a listing  of  the  code  which  includes 
comment  cards.  This  section  covers  aspects  relating  to  the  logic, 
structure,  and  operating  requirements  of  the  entire  program. 

5.1  PROGRAM  LOGIC 

The  program  logic  is  described  by  its  natural  division  into  initial- 
ization and  reinitialization,  data  input,  preprocess  and  error  checking, 
computation,  and  output.  Table  5 -A  shows  how  the  FAST  routines 
divide  into  the  various  functions.  In  addition,  Figure  5-1  outlines  the 
logic  flow  or  calling  sequences  and  identifies  which  subsection  (of  5.  1) 
describes  each  calling  sequence. 

As  Figure  5-1  shows,  the  program  begins  with  the  routine  called 
MAIN,  which  calls  subroutine  NLREAD.  During  execution,  NLREAD 
drives  the  program  logic  by  initiating  calls  to  various  routines  or 
sequences  of  routines.  When  a STOP$  card  is  encountered  at  the  end 
of  the  data,  NLREAD  returns  to  MAIN  where  the  execution  is  terminated. 

The  logic  in  NLREAD  is  primarily  controlled  by  the  names,  like 
STOP$,  that  are  punched  on  the  data  cards.  These  names  as  well  as 
data  names  are  described  in  Section  3.  0 and  can  be  seen  in  the  listing 
of  MAIN.  MAIN  stores  the  names  in  common  block/NLIST/  which  is 
passed  to  subroutine  NLREAD. 

The  following  description  of  the  program  logic  assumes  the 
reader  will  refer  to  or  be  familiar  with  the  input  (particularly  the  input 
names),  the  listing  of  NLREAD,  Table  5-A  and  Figure  5-1. 

5.1.1  Initialization  and  Reinitialization 

The  program  execution  begins  in  MAIN,  which  initializes  tables 
and  variables  that  need  to  be  set  only  once  for  a computer  run  (MAIN 
serves  in  place  of  a block  data  routine  since  this  feature  is  not  avail- 
able on  all  computers).  The  rest  of  the  variables  are  initialized  by 
subroutine  RESET,  which  is  automatically  called,  by  NLREAD,  before 
the  program  reads  the  first  job. 
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INITIALIZATION  AND  REINITIALIZATION 


igure  5-1.  Hierarchy  of  Subroutine  Call 
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Table  5-A.  Abstracts  of  FAST  Routines 


* Initialization  and  Reinitialization 


MAIN 


RESET 


Data  Entry 


NLREAD 


NLGET 


NLNAME 


NLNUMB 


NLREAL 


NLINTG 


NLVECR 


NLVECI 


COMPNT 


Interfaces  program  with  different  computer 
systems,  dimensions  program  vectors  and 
matrices,  and  serves  as  a block  data  routine 
to  initialize  tables  and  variables  that  are  set 
only  once  for  a computer  run. 


Initializes  those  variables  that  must  be  set  before 
each  job  in  a computer  run. 


Controls  the  program  flow  by  calling  the  appro- 
priate routines  for  each  input  name  read  by 
NLNAME. 


Reads  and  prints  80  column  data  cards.  NLGET 
also  provides  the  card  character  recode  option. 


Identifies  input  names  (read  from  data  cards)  and 
converts  subscripts  on  the  cards  to  integer 
numbers. 


Converts  card  numbers  to  signed  integer,  fraction 
and  exponent.  NLNUMB  also  controls  how  many 
times  the  number  is  to  be  repeated,  identifies 
when  NLREAL  should  take  log  base  10  of  the  num- 
ber, and  whether  the  next  card  field  contains  a 
number  or  input  name. 


Forms  and  stores  real  (floating  point)  numbers 
read  by  NLNUMB. 


Forms  and  stores  integer  numbers  read  by 
NLNUMB. 


Repeatedly  calls  NLREAD  to  enter  a vector  of 
real  numbers.  NLVECR  indexes  and  counts  the 
vector  elements. 


Repeatedly  calls  NLINTG  to  enter  a vector  of 
integer  numbers.  NLVECI  indexes  and  counts 
the  vector  elements. 


Interprets  and  stores  component  card  information. 
COMPNT  is  also  used  to  delete  a component  for 
case  stacking. 
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Table  5-A.  Abstracts  of  FAST  Routines  (Continued) 


Data  Entry  (Continued) 

BOOLE  Stores  a complete  network  equation  from  the 

equal  (=)  symbol  to  the  dollar  ($)  symbol;  trans- 
lates the  component,  subsystem,  and  system  names 
to  addresses;  and  calls  subroutine  LOGIC.  BOOLE 
is  also  used  to  delete  equations  for  case  stacking. 

LOGIC  Translates  the  network  equation  into  a symbolic 

table  of  instructions  for  finding  the  probability  of 
failure. 

Preprocess  and  Error  Checking 

RTOP  Converts  ranges  to  equivalent  pressures. 

PTOR  Converts  pressures  to  equivalent  ranges. 

SELECT  Selects  from  the  data  base  those  environments 

required  for  running  a case  and  scales  the  ULN  to 
yield  and  height  of  burst. 

NETCK  Checks  the  network  table  of  instructions  for 

possible  errors. 


Computation 

PROG  Performs  the  Monte  Carlo  assessment  of  compo- 

nent, subsystem,  and  system  probability  of  failure. 

RNVD  Decomposes  the  covariance  matrices. 

CAN  Computes  the  sum  of  environments  as  part  of  the 

transfer  function.  This  routine  will  sum  log 
base  10  environments  by  taking  antilogs  (this  rou- 
tine cannot  combine  log  and  non-log  environments). 

USER  As  part  of  the  input,  the  user  can  provide  this 

routine  for  azimuth  sensitivity. 

TABLE  If  the  azimuth  data  is  in  tabular  form  USER  can 

call  TABLE  - no  other  FAST  routines  call  TABLE. 
TABLE  performs  a table  look-up  with  linear 
inte  rpolation. 

[included  in  the  computation  section  are  two  random  number 
generators  - one  for  uniformly  distributed  numbers  and  one 
for  normally  distributed  numbers.  These  routines  are  not 
identified  by  name  because  the  names  will  be  different  on 
other  computer  systems.  Sec  comment  cards  in  subroutine 
PROG.  ] 
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Table  5-A.  Abstracts  of  FAST  Routines  (Continued) 


Output 

PROMED 

Prints  the  histogram  and  processes  the  histogram 
to  find  the  median  and  other  confidence  levels  for 
printing.  The  medians  are  optionally  saved  for 
DOPLOT. 

PROBL 

Used  by  both  PROG  and  PROMED  to  integrate 
histogram  data  to  find  median  and  confidence 
levels  for  convergence  tests  and  printing. 

RORD 

Reorders  the  components,  subsystems,  and  sys- 
tems by  hardness  of  the  median  values. 

DOPLOT 

Forms  plot  borders  and  f'^ntrols  plotting 
operation. 

SPLOT 

Determines  plot  symbol  and  position  for  each 
point  on  plot  page. 

PRPLOT 

Prints  lines  with  two  plots  per  page. 

A run  may  contain  more  than  one  job  where  each  job  is  like  a 
separate  computer  run.  A RESET$  card  identifies  the  beginning  of  a 
new  job  and  causes  NLREAD  to  again  call  subroutine  RESET. 


There  is  one  important  difference  between  separate  computer 
runs  and  stacked  jobs  - the  starting  values  (NUMUR  and  NUMNR)  for 
the  random  number  generators  are  not  reset  for  each  job.  This  means 
that  successive  jobs  in  a run  use  different  sequences  of  random  num- 
bers, whereas  separate  runs  start  wifh  the  same  sequence  each  time. 
However,  the  RANUN  data  card  n ay  (on  some  computer  systems)  be 
used  to  set  or  reset  NUMUR  and  NUMNR. 


There  are  some  other  variables  that  are  not  reset  in  job  stacking. 
These  variables,  which  are  affected  by  CONF1,  PLOTON  and  RECOPE 
data  cards,  can  be  found  in  the  listing  for  MAIN  under  the  comments 
"SET  DEFAULT  VALUES"  and  "INITIALIZE  READING  OF  CARDS  BY 
NLGET".  Generally,  all  jobs  in  a run  will  use  the  same  conditions  so 
not  resetting  these  variables  is  not  of  concern. 
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Job  stacking  should  not  be  confused  with  case  stacking.  The 
RESET$  card  purges  all  previous  information,  but  case  stacking  only 
replaces  or  removes  selected  information,  which  is  done  in  the  data 
entry  routines. 


5.1.2  Data  Entry 

Once  the  program  has  been  initialized  by  MAIN  and  RESET, 
NLREAD  calls  NLGET  to  read  the  first  data  card.  NLGET  is  the  only 
routine  that  reads  cards,  and  is  the  only  routine  that  recognizes  and 
processes  the  RECODE  card  (the  RECODE  card  is  not  returned  as  a 
data  card).  Before  an  80-character  data  card  is  printed  and  returned 
in  the  80-word  vector  called  ICARD,  the  characters  are  recoded  if  the 
recode  option  is  in  effect. 


The  call  to  NLGET  is  within  the  first  of  two  major  loops  in 
NLREAD  and  is  executed  after  a RUN$  or  PLOT$  card,  where  it  is 
necessary  to  print  other  information  before  the  next  data  card  is  read 
and  printed.  If  in  the  future  it  is  necessary  to  read  fixed  format  data, 
a new  name$  could  be  introduced;  then  following  the  name$  would  be  the 
fixed  format  cards  which  would  be  processed  before  a new  call  to 
NLGET. 


The  second  major  loop  in  NLREAD  begins  with  a call  to  NLNAME. 
Since  the  cards  are  punched  in  free-field  format,  the  data  are  treated 
as  a sequence  of  logical  records  where  each  record  begins  with  an  iden- 
tifying input  name.  NLNAME  is  used  to  translate  each  name  into  an 
integer  code  NAMEIX,  which  is  used  by  NLREAD  to  branch  to  a section 
of  code  to  process  the  logical  record. 


The  data  entry  routines  process  logical  records  that  contain  a list 
of  information  to  be  read  and  stored.  The  forma  for  two  of  the  logical 
records  is  sufficiently  complex  to  require  special  routines  - NLREAD 
calls  COMPNT  to  process  the  component  cards  and  BOOLE  to  process 
the  network  cards.  The  format  for  the  rest  of  the  data  is  embedded  in 
NLREAD  which  makes  use  of  tb  » generalized  routines  NLVECI, 
NLVECR,  NLINTG  and  NLREAL  (see  Table  5-A). 
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Subroutine  NLNUMB  is  used  to  read  numeric  fields  on  the  cards. 

Every  field  is  assumed  to  have  the  following  format: 

IT 

r * ± i . f ^ ± eT 

NLNUMB  stores  the  various  parts  of  this  format  in  common  block/ 

CARDAT /where  NLINTG  and  NLREAL  form  the  parts  into  a number. 

r*  The  * identifies  and  separates  r from  the  rest  of  the  field. 

The  r is  stored  as  IREP  and  is  the  number  (integer)  of 
times  NLNUMB  will  be  called  before  NLNUMB  formats 
the  next  field  (3  * 4$  is  the  same  as  4,  4,  4 $). 

! 

±i  . f The  positive  (the  + is  ignored  and  need  not  be  punched)  or 
negative  number  is  stored  in  two  parts  in  INTPRT  = (±  i) 
and  FRCPRT  - (±  • f).  The  number  of  digits  in  i and  f 
are  added  together  and  stored  in  NULCNT.  If  NULCNT  is 
non-zero,  then  NLINTG  or  NLREAD  will  form  and  store 
the  new  number  (in  place  of  a number  set  by  MAIN  or 
RESET  or  entered  on  a previous  card).  If  NULCNT  = 0 
then  the  number  which  was  read  in  a previous  case  or  on 
some  other  cards  is  not  replaced. 

i 

E 

^ ± e The  E and  L separates  the  exponent  from  the  number  and 
the  positive  or  negative  integer  value  of  e is  stored  in 
IEXPRT.  Again,  the  + does  not  need  to  be  punched.  If  X 

4 0 

represents  the  number,  then  E means  X*  10  where  ILOG 
will  be  set  to  zero.  The  L means  to  take  the  log  base  10 
of  X*  1 0±G  or  ±e  + log  10  (X),  where  ILOG  - 1. 

T T terminates  or  separates  one  field  from  another.  T is 

interpreted  by  NLNUMB  and  represented  by  the  integer 
code  NEXT.  Since  NEXT  is  a complicated  parameter  to 
define,  the  conditions  are  illustrated  in  Table  5-B. 


70 





able  5-B.  Illustration  of  Variable  NEXT 
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R1=(NEXT=4  - non-numeric  equal  read  by  COMPNT) 
3$  (NEXT  =j  - number  followed  by  dollar  symbol) 
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There  is  an  execution  inhibit  flag  (INHIBT)  in  common  block 
/CARDAT/  that  is  incremented  (for  each  error)  if  the  program  finds  an 
error  in  the  data.  This  flag  prevents  the  program  from  running  the  job 
and  machine  time  is  not  wasted  by  trying  to  perform  the  Monte  Carlo 
assessment.  If  a RESET$  card  is  found,  subroutine  RESET  turns  the 
flag  OFF  (INHIBT  = 0),  because  the  next  job  will  have  a completely 
new  set  of  information.  In  order  to  continue  scanning  the  data  NLNUMB 
is  used  to  skip  numeric  fields  until  the  next  name  is  found,  which  starts 
the  next  logical  record. 


5.  1.  3 Preprocess  and  Error  Checking 

The  RUN$  card  indicates  that  all  of  the  data  has  been  read  for  a 
case.  At  this  point  NLREAD  checks  the  inhibit  flag.  If  it  is  non-zero 
(which  means  there  is  an  error)  NLREAD  transfers  to  the  major  loop 
(in  NLREAD)  that  reads  the  next  data  card  and  the  program  remains  in 
the  data  entry  procedure.  If  the  inhibit  flag  is  zero  then  NLREAD  calls 
the  preprocess  and  error  checking  routines  SELECT  and  NETCK  before 
calling  PROG  to  perform  the  Monte  Carlo  assessment.  Since  SELECT 
and  NETCK  can  also  set  the  inhibit  flag,  the  flag  is  again  checked 
before  calling  PROG  and  the  same  transfer  takes  place  as  was  just 
described. 


FAST  has  an  environment  "data  bank"  capability.  This  means 
that  the  user  may  enter  more  environments  than  are  required  for  the 
case.  Since  the  component  cards  identify  which  environments  are 
required,  subroutine  COMPNT  maintains  a list  of  the  environments  in 
the  partitioned  vector  called  MODE  by  adding  a one  to  an  element  of 
MODE  each  time  the  environment  is  requested.  If  a component  is 
deleted  or  replaced,  COMPNT  subtracts  one  from  the  appropriate 
elements  of  MODE. 


Subroutine  SELECT  scans  the  MODE  vector  to  determine  which 
environments  are  required  and  then  condenses  the  environment  informa- 
tion (particularly  the  covariance  matrix  and  correlation  matrix). 
SELECT  also  scales  the  input  ULN  from  the  nominal  yield  and  height 
of  burst  to  the  yield  and  height  of  burst  of  interest. 
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Both  SELECT  and  NLREAD  call  the  range  to  pressure  (RTOP) 
and  pressure  to  range  (PTOR)  routines.  This  allows  the  user  to  enter 
either  ranges  or  pressures  and  the  program  will  compute  the  equivalent 
values,  since  both  are  required  by  the  program. 

NETCK  scans  the  Boolean  networks  for  errors.  Several  error 
messages  can  be  printed  which  are  described  in  Section  4.  1 and  can  be 
found  in  the  NETCK  listing. 

5.1.4  Computation 

If  no  errors  were  found  in  the  data,  a RUN$  card  causes  NLREAD 
to  call  subroutine  PROG  to  perform  the  Monte  Carlo  assessment. 
Although  PROG  is  called,  the  program  returns  to  NLREAD  if  RNVD 
finds  that  a covariance  matrix  is  not  positive  semidefinite  (see 
Section  2.  2.  5.  2). 

The  reason  RNVD  is  called  from  PROG  rather  than  SELECT,  is 
that  the  K-factor  uncertainty  for  log  base  10  environments  is  dependent 
on  the  scaled  mean.  Thus,  the  uncertainty  covariance  must  be  con- 
structed from  a new  K-factor  each  time  the  ULN  is  called.  Log  Base  10 
environments  are  not  commonly  used,  so  this  process  is  not  discussed 
in  the  text  although  a programmer  could  determine  what  is  done  from 
the  PROG  listing. 

PROG  performs  the  Monte  Carlo  assessment  of  component,  sub- 
system and  system  survivability.  This  process  is  described  by  com- 
ment cards  in  the  PROG  listing  showing  the  major  loops  and  the  purpose 
of  each  section  of  code.  Because  the  program  is  designed  to  simplify 
redimensioning  variables,  which  is  discussed  in  Section  5.2.  1,  the  way 
information  is  stored  and  indexed  is  somewhat  complicated. 

It  is  beyond  the  scope  of  this  manual  to  trace  each  variable  through 

y 

the  program  and  show  how  it  is  used,  but  the  following  information  may 
be  useful.  Section  5.  2.  2 describes  the  common  blocks  as  to  how  most 
of  the  data  is  stored  and  what  information  the  blocks  contain.  The  envi- 
ronment data  is  read  by  NLREAD,  rearranged  in  SELECT  and  used  by 
PROG.  The  component  cards,  including  fragility  curves  and  transfer 
functions,  are  read  by  COMPNT  and  used  by  PROG.  Since  component 
cards  may  be  deleted  or  replaced,  they  are  not  necessarily  stored  in 


the  same  order  as  they  are  read,  but  can  be  located  by  indirect  indexing 
using  common  block  /CC01/.  The  same  is  true  of  the  network  data, 
which  can  be  found  using  common  block  /CNC2/.  The  network  is  read 
by  BOOLE  and  the  translated  instructions  are  stored  by  LOGIC  in  com- 
mon block  /CN01/  which  is  used  by  PROG. 

5.1.5  Output 

Subroutine  PROMED,  which  is  called  by  PROG,  prints  the  histo- 
gram and  confidence  matrix.  These  results  are  printed  for  each  range/ 
pressure  that  is  read  from  either  the  R or  OP  data  cards.  If  the  plot 
option  (PLOTON)  has  been  requested,  PROMED  stores  the  50%  values  in 
the  matrix  PROB  of  common  block  /CT07/,  which  is  used  by  the  printer 
plot  routines. 

N1  READ  calls  DOPLOT  when  a PLOT$  card  is  encountered  in  the 
datu.  Since  cases  may  have  to  be  stacked  to  change  the  scaling  param- 
eters as  the  value  of  pressu-e  increases,  the  PLOT$  card  may  appear 
after  the  last  case  and  causes  all  of  the  results  to  be  plotted  as  a single 
plot.  This  is  accomplished  since  common  block  /CT07/  can  contain 
results  from  several  RUN$  cards. 

In  order  to  have  a machine  independent  plot  capability,  the  printer 
is  used  to  plot  the  data.  The  plots  are  sufficiently  accurate  unless  the 
available  printer  produces  wavy  lines.  Currently  only  the  50%  curves 
for  subsystems  and  systems  are  sorted  by  PROMED  and  plotted  by 
DOPLOT,  SPLOT  and  PRPLOT. 

5.2  PROGRAM  STRUCTURE 

The  program  is  structured  to  provide  maximum  problem  solving 
flexibility  to  the  user  while  maintaining  ease  of  usage.  On  many  compu- 
ters it  is  cost  effective,  with  a long  running  Monte  Carlo  program,  to 
use  a minimum  amount  of  core  memory.  For  this  reason,  care  has 
been  taken  to  simplify  the  redimensioning  of  variables  when  unusually 
large  or  complex  problems  must  be  solved.  FAST  is  used  on  various 
computers  so  a special  emphasis  has  been  placed  on  making  the  program 
as  machine  independent  as  possible.  These  features,  along  with  com- 
mon block  organization,  are  discussed  in  the  following  sections. 
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5.2.1  Variable  Redimensioning 


Most  variables  are  completely  dimensioned  only  in  the  main  driver, 
and  are  dimensioned  as  dummy,  one-word  vectors  in  all  other  routines. 
Dimensions  depend  on  the  size  of  the  problem  to  be  solved,  such  as  the 
number  of  components,  the  number  and  complexity  of  network  equations, 
or  the  number  of  environments.  In  redimensioning  the  program,  it  is 
necessary  to  make  sure  that  the  maximum  problem  size  is  properly  set 
by  MAIN  and  that  interrelated  as  well  as  directly  related  dimensions  are 
changed.  Table  5-C  and  the  listing  for  MAIN  better  illustrate  what  needs 
to  be  done,  and  if  a mistake  is  made  there  are  some  error  messages  that 
may  be  helpful  (see  Section  4.  2). 

The  technique  used  to  simplify  redimensioning  makes  use  of 
properties  of  common  blocks.  As  the  program  is  loaded,  the  loader 
allocates  core  by  the  maximum  dimension  of  the  common  block  in  MAIN. 
In  the  subroutines  the  blocks  appear  shorter  but  the  starting  address 
and  amount  of  available  core  are  the  same.  The  dimensioned  variable 
is  always  the  last  parameter  of  the  block.  It  should  be  noted  that  the 
variables  IPAD1,  IPAD2  and  IPAD3  are  not  used  by  the  program,  but 
are  used  to  force  real  variables  to  even  word  boundaries  for  the  IBM  360. 

5.2.2  Common  Blocks 

The  letters  in  the  common  block  names  help  to  identify  what  kind 
of  information  is  stored  in  the  block  {see  Table  5-C).  Table  5-D  cross 
references  common  blocks  to  subroutines.  As  mentioned  in  Sec- 
tion 5.2.  1,  Table  5-C  indicates  how  the  blocks  are  dimensioned. 

The  environment  common  blocks  are  more  complex  than  the  others, 
because  the  program  has  an  environment  data  bank  capability.  The  read 
routines  store  all  the  possible  environments  in  an  alternate  region. 

Then,  depending  on  which  environments  are  called  for  by  the  component 
cards,  subroutine  SELECT  copies  those  environments  to  the  primary 
region  prior  to  the  Monte  Carlo  iterations.  The  data  bank  capability 
reduces  the  chance  that  an  environment  will  be  missing,  and  at  the  same 
time,  it  keeps  the  program  from  performing  unnecessary  computations. 
The  MODE  vector,  which  is  partitioned  into  four  parts,  contains  the 
selection  information  and  is  used  by  subroutine  SELECT. 
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Table  5-C.  Common  Blocks  - Contents  and  Dimensions 
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FRCPRT  Fractional  part  of  the  number. 


Table  5-C.  Common  Blocks  - Contents  ar-1  Dimensions  (Continued) 
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enerally  only  50%  (.5)  confidence  is  tested 
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5.2.3  Computer  Conversion 

Program  conversion  to  other  computers  is  facilitated  by  requiring 
changes  only  to  MAIN  and  subroutine  PROG.  All  other  routines  are 
believed  to  be  machine  independent. 

MAIN  defines  the  FORTRAN  logical  input/output  units,  which  are 
easy  to  change.  MAIN  also  contains  all  data  statements  because  the 
format  may  be  different  on  some  computers  (e.g.  , the  CDC  3800). 


Since  routines  to  generate  random  numbers  (see  Section  5.  3)  are 
machine  dependent,  the  calls  to  the  routines  in  subroutine  PROG  need 
to  be  changed.  If  the  generators  require  starting  values,  the  variables 
NUMUR  and  NUMNR  can  be  set  in  the  main  driver  or  by  placing  a 
RANUN  card  in  the  data. 


The  program  is  structured  to  run  on  ar  IBM  360  or  similar  IBM 
computer,  and  even  provides  for  the  possibility  that  a data  deck  may 
contain  cards  punched  on  both  an  026  and  029  keypunch  (see  the  descrip- 
tions for  the  RECODE  data  card). 


5.  3 OPERATING  REQUIREMENTS 

The  FAST  program  requires  a card  input  unit  and  a printer  output 


init.  No  other  input/output  devices  are  currently  required. 


On  the  CDC  computers,  it  is  necessary  to  load  MAIN  first  since 
this  is  the  only  routine  that  dimensions  common  block  variables  - this 
avoids  truncated  common  blocks. 


The  program  requires  two  random  number  generators,  one  that 
produces  numbers  uniformly  distributed  between  0 and  1,  and  one  that 
produces  numbers  from  a normal  distribution  with  mean  of  zero  and 
variance  of  unity. 


The  core  storage  requirements  depend  upon  the  efficiency  of  the 
loader  and  compiler  (the  TRW  FUN  optimizing  compiler  lias  been  used 
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successfully).  Although  FAST  is  used  on  various  computers,  the  only 
current  core  storage  information  is  for  TRW  and  Lawrence  Berkeley 
Laboratory. 


LOAD 


EXECUTE 


CDC  6600 
TRW 

CDC  7600 
LBL 

36  K8 

73  Kg 

66  K8 

63  K8 
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6.0  SAMPLE  PROBLEM 


A sample  problem  has  been  included  in  this  manual  to  provide  a 
check  case  for  comparing  output,  to  accelerate  the  user's  acquisition  of 
insight  to  FAST  techniques,  to  develop  intuition  as  to  what  constitutes  a 
reasonable  output  from  FAST  given  a reasonable  set  of  inputs,  and  to 
describe  the  development  of  FAST  inputs  for  a "real"  physical 
configuration. 


To  accomplish  these  purposes,  FAST  inputs  for  performing  surviv- 
ability studies  for  a military  system  subject  to  a nuclear  threat  have  been 
developed.  The  sample  problem  considers  only  effects  due  to  the  air- 
blast,  ground  shock,  and  crater  debris  environments,  however,  other 
nuclear  weapon  environments  such  as  EMP,  nuclear  and  thermal  radia- 
tion, etc.,  could  be  considered  within  the  FAST  methodology.  The  sample 
problem  was  designed  to  exhibit  important  aspects  of  modeling  the  system 
and  the  response  phenomenology  in  FAST  and  considers  a hypothetical  sys- 
tem and  generic  site  conditions.  The  different  aspects  of  developing  the 
FAST  inputs  are  described  in  varying  degrees  of  detail,  depending  upon 
their  significance  in  demonstrating  the  code  capability  and  practical 
modeling  considerations.  In  particular,  the  development  of  nuclear  blast 
and  shock  environment  inputs  has  been  emphasized,  with  highlights  on  the 
shock  spectral  approach  for  evaluating  ground  shock  induced  damage. 


The  hypothetical  system  for  the  sample  problem  consists  of  an 
ICBM  control  facility  shown  in  Figure  6-1.  A number  of  these  systems, 
built  almost  identically,  are  to  be  deployed  at  sites  where  geology  is 


soil  with  a water  table  at  least  200  feet  below  the  surface  and  bedrock 


below  the  water  table.  The  design  hardness  is  600  psi  for  a five  megaton 
(5  MT)  surface  burst.  The  system  is  in  an  early  stage  of  development, 
and  the  survivability  of  a baseline  system  is  to  be  evaluated.  This  study 
will  consider  the  system  survivability  as  a function  of  overpressure  for 
yields  of  0.  5 and  5 MT . and  both  surface  contact  and  near-surface  height 


of  burst  detonations. 


To  perform  its  mission,  the  control  facility  must  be  able  to  receive 
and  transmit  messages  by  using  either  the  communications  anter.na  or  the 
hardened  cable  link  (HCL).  The  command  and  control  consoles,  which  are 
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essential  to  mission  performance,  are  powered  by  batteries  and  cooled  by 

* 

air  conditioning.  A sturdy  rectangular  structure  consisting  of  two  nearly 
cubical  rooms  protects  the  equipment.  In  normal  operation,  the  control 
facility  circulates  air  from  the  outside.  When  attacked,  the  blast  valve 
closes  the  air  entrainment  system,  thus  protecting  the  control  room. 

Table  6-A  shows  a detailed  hardness  matrix  for  the  system.  The 
matrix,  which  relates  the  various  functional  subsystems  and  components 
with  the  environments  to  which  they  are  most  sensitive,  constitutes  a 
compact  means  of  controlling  FAST  inputs.  The  numbers  in  the  environ- 
ment column  refer  to  the  corresponding  parameters  (listed  in  i able  6-C). 

6.  1 ENVIRONMENTS 

This  section  deals  with  the  derivation  of  hostile  environment  inputs 
for  the  FAST  code  sample  problem.  For  this  hypothetical  system,  a 
design  threat  of  a surface  burst  of  5 MT  yield  at  a range  corresponding  to 
600  psi  has  been  chosen  as  the  baseline  condition  from  which  weapon  effects 
predictions  are  to  be  scaled.  Because  of  nonlinearities  in  scaling  over 
large  distances,  the  baseline  threat  conditions  should  be  as  near  to  the 
middle  of  the  range  of  interest  as  possible. 

Ground  shock  environment  prediction  equations,  shown  in  Table  6-B, 
were  developed  in  coordination  with,  and  under  the  direction  of  DNA/RDA, 
and  provide  the  basis  for  derivation  of  the  environments.  These  equations 
consider  the  peak  free  field  air  induced  vertical  and  horizontal  displace- 
ments, velocities,  and  accelerations  and  the  crater  and  related  effects 
(CARE)  displacement  and  velocity  (where  the  vertical  and  horizontal  CARE 
components  are  taken  to  be  equal).  In  the  equations,  L is  the  minimum  of 
the  depths  to  bedrock  and  water  table  (L  = 300  ft),  z is  the  free  field  depth 
of  interest,  t^  is  the  seismic  arrival  time  from  the  surface  to  the  depth 
of  interest,  p is  the  media  density,  U is  the  airblast  shock  velocity,  and 
V is  the  crater  volume  (v  = 50  x 10^  ft^  | [ ) • The  free  field  vertical 

shock  spectra  for  the  25-foot  depth  and  the  free  field  horizontal  shock 
spectra  for  the  15-foot  depth  (for  the  design  threat)  are  shown  in  Fig- 
ure 6-2.  The  various  aspects  of  this  spectra  were  designed  partly  to  show 
the  detail  to  which  the  spectra  can  be  modeled  and  does  not  necessarily 
reflect  the  detail  required  for  a given  problem.  The  lower  frequency 
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Table  6-A.  Sample  Problem  Hardness  Matrix 


Table  6-B.  Ground  Shock  Environment  Equations  - Surface  Burst 


response  of  the  structure  would  be  governed  by  these  spectra.  The 
median  values  for  the  baseline  threat  condition  and  the  scaling  laws  are 
shown  in  Table  6-C.  In  that  table,  the  parameters  fall  into  four 
categories: 

a.  Free  field  peak  air-inducec]  ground  shock  displacements, 
velocities  and  accelerations 


i 

I 


| 

I 


! 
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b.  Derived  air- induced  ground  shock  parameters  which 
primarily  specify  line  segments  on  the  shock  spectra 
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Table  6-C.  Median  Environment  at  5 MT,  600  psi  plus  Scaling  Laws 

Scaling 


Environment 

Parameter 


Median 

Value* 


E~-R  P W° 


Static  Footing  Load 
E(17)  dCARE 


E(18) 


CARE 

E(19)  Debris  Depth 


27  in. 
29  ips 
12  in. 


■ 3.0  1.33 

■ 2.0  2/3 


-3.86  1.6 


*5  MT,  600  psi  - Surface  Burst 


-n^w0 


E-p’w 


P 

a 

P 

a 

E(l) 

d 

z 

16  in. 

2/3 

? /6 

E(2) 

V 

z 

3 80  ips 

1.0 

1/12 

E (3 ) 

az 

60  g's 

1.0 

1/12 

E(4) 

V 

z max 

570  ips 

1.0 

1/  12 

E(5) 

20  d 2/v 
z z 

14  in.  sec 

1/3 

1/4 

E ( 6) 

dH 

4 in. 

2/3 

1/6 

E ( 7 ) 

VH 

45  ips 

1/2 

1/12 

E(8) 

Vrr 

H max; 

67  ips 

1/2 

1/12 

E(9) 

1.5  vH2/dH 

1.9  g*  s 

1/3 

0 

E(10) 

'‘'rel 

2400  lb/in. 

1.8 

1/12 

E(1 1) 

Dynamic 

Pressure 

1400  psi 

1.37 

0 

E(12) 

agTR(A/C) 

10C  g’s 

1.0 

1/12 

E ( 1 3 ) 

aSTR^SSS^ 

200  g's 

1.0 

1/25 

E(14) 

Equivalent 
Static  Roof  Load 

600  psi 

1.0 

1/25 

E(  1 5) 

Equivalent 
Static  Walls  Load 

600  psi 

. 1-0 

1/12 

E(16) 

Equivalent 

600  psi 

1.0 

1/6 

H 
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c.  Correlated  air-induced  effects  on  the  control  facility 
structure 

d.  Cratering  and  Related  Effects  (CARE). 

Table  6-C  shows  that  air-induced  environments  scale  with  over- 
pressure and  yield,  while  crater-induced  effects  scale  with  range  and 
yield.  For  surface  bursts  on  a given  geology,  CARE  effects  scale  with 
crater  volume  and,  thusly,  with  yield.  For  heights  of  burst,  the  crater 
volume  reduction  affects  the  CARE  environments  while  changes  in  the 
overpressure/range  law  at  different  heights  of  burst  are  used  to  estimate 
changes  in  Pressure  and  Related  Environments  (PARE). 

The  list  of  environments  in  Table  6-C  are  designed  to  model  the 

critical  system  response  phenomenology  in  a statistically  meaningful 

manner.  Heading  the  list  of  weapons  effects  are  iree  field  vertical  peak 

displacement,  d , velocity,  v , and  acceleration,  a , at  the  25-foot  depth 
z z z 

below  the  surface.  These  parameters  characterize  the  rigid-body  vertical 
motion  of  the  control  facility.  Next  comes  V^,  the  peak  vertical  shock- 
spectral  velocity  (shown  in  Figure  6-2). 

Along  the  displacement  amplification  segment  of  the  spectrum  the 

displacement-to-frequency  ratio  (D/f)  is  constant.  The  pseudo- velocity 

2 

(V)  must  then  be  of  the  form,  V = pf  , where  the  proportionality  factor  (p) 
is  evaluated  from  the  pseudo- velocity  at  the  break  point  frequency  (see 
Figure  6-2) 


f = To  cl  CPS  an*ts) 


(6-1) 


The  break  point  pseudo-velocity  is  equal  to  (tt/10)  v and  therefore  the 
displacement  amplification  segment  can  be  described  by 


V = 40Tri 


(6-2) 


The  environment  parameter 


20  d /v 
z z 


(6-3) 
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specifies  the  displacement  amplification  segment  of  the  shock  spectrum  in 
this  analysis.  The  horizontal  peak  displacement,  d. and  vertical  peak 
velocity,  v^,  in  the  free  field  at  the  15-foot  depth  characterize  the  hori- 
zontal rigid-body  response  of  the  control  facility.  The  horizontal  peak 

2 

shock  spectral  velocity  is  V^.  The  quantity  1.  5 ilmes  the 

constant  acceleration  line  connecting  the  two  constant  pseudo-velocity 
lines.  The  quantity  o"dre^  the  critical  parameter  for  the  cable  pene- 
tration response.  It  is  the  product  of  the  backfill  stress  times  the  back- 
fill motion,  relative  to  *he  control  facility.  The  pressure,  Q,  is  the  peak 
dynamic  pressure  associated  with  the  passage  of  an  overpressure  shock. 
The  structural  acceleration  environments  are  given  at  the  ceiling  attach- 
ment points  of  the  shock  isolated  floor,  ag^.^  (SSS),  and  at  the  location  of 


the  air  conditioner,  ag^,^  (A/C).  Note  that  these  environments  relate 


more  to  the  structural  configuration  than  to  the  geological  site  condition. 
The  final  two  classes  of  environments  relate  to  structural  collapse  loads 
and  to  CARE  environments.  The  CARE  environments,  displacement, 
velocity  and  debris,  are  self-explanatory.  However,  the  structural  load 
environments  require  further  explanation. 


The  three  structural  subsystems  (roof,  walls,  and  footing)  each 
have  fragilities  which  are  functions  of  parameters  closely  related  to  over- 
pressure. It  is  important  to  consider  correlation  in  tKe  analysis  of  these 
subsystems;  for  example,  soft  foundation  materials  tend  to  cause  a rela- 
tively large  foundation  loading  and  cushion  the  roof  loading,  while  for  the 
case  of  stiff  foundation  materials,  the  reverse  is  true.  The  key  construc- 
tion variables  relating  to  the  collapse  of  structural  members  include  con- 
crete strength,  and  workmanship  in  placing  steel,  concrete,  and  backfill. 
These  construction  variables  also  introduce  correlation  in  the  structural 
response.  The  FAST  code  does  not  permit  correlation  of  the  fragilities; 
however,  the  proper  correlations  between  the  structural  subsystems,  due 
to  construction  techniques  and  structural  mechanics,  can  be  obtained 
through  correlation  of  the  structural  loads. 


The  environments  giving  equivalent  static  roof,  wall  and  footing 
loads  take  into  account  the  dynamic  elasto-plastic  responses  in  the  failure 
modes  of  these  structural  components.  The  scaling  with  yield  indicates 
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that  the  roof  and  walls  which  respond  in  shorter  times  than  the  footing  are 
less  yield-dependent  than  the  footing,  i.  e.  , the  footing  load  must  be  sus- 
tained for  a longer  time  in  order  to  produce  collapse.  Therefore,  the 
footing  response  is  more  yield-dependent.  The  impulse,  I = Jq  P(t)dt,  up 
to  the  critical  time  for  structural  collapse,  t*,  is  the  key  variable  and  is 
accounted  for  in  the  yield  scaling  factor. 

Now  that  all  the  environment  parameters  have  been  addressed,  it  is 
worthwhile  to  further  discuss  the  advantages  of  the  shock  spectrum 
approach  in  hardness  analysis.  The  use  of  shock  spectra  greatly  simplifies 
the  prediction  of  the  damage  potential  of  various  ground  shock  time  his- 
tories. The  shock  spectrum  of  a shock  waveform  is  the  peak  response  to 
the  waveform  of  a simple  oscillator  as  a function  of  the  oscillators'  fre- 
quency. Figure  6-2  shows  two  shock  spectra.  The  usefulness  of  shock 
spectra  in  damage  prediction  stems  from: 

a.  The  peak  response,  for  many  structural/rnechanical 
systems,  can  be  determined  directly  from  the  shock 
spectra. 

b.  Shock  waveforms  of  the  same  class  often  will  have  much 
the  same  shock  spectra  even  though  the  actual  time 
histories  differ  considerably. 

Given  the  orientation  of  a line  on  shock  spectra  paper,  only  one 
parameter  is  needed  to  specify  the  line.  For  example,  the  solid  line  seg- 
ments below  0.2  cps  in  Figure  6-2  are  lines  of  constant  maximum  relative 
displacement,  D;  the  lines  above  100  cps  are  lines  of  constant  maximum 
acceleration  response,  A;  constant  pseudo- velocity,  V,  lines  are  horizon- 
tal. For  any  point  on  this  type  of  paper,  D,  V and  A,  as  well  as  frequency, 
are  specified.  Over  the  years  engineers  have,  for  certain  classes  of 
pulses,  used  simple  but  accurate  rules  for  estimating  shock- spectral  dis- 
placement, D = d,  shock- spectral  peak  pseudo-velocity,  V = 1.5v,  and 
shock-spectral  acceleration,  A = 2a,  where  capital  letters  are  spectrum 
parameters  and  small  letters  are  free  < old  values  of  peak  displacement, 
velocity  and  acceleration.  In  general,  the  shock  response  spectrum  values 
are  amplified  over  peak  free  field  values.  More  recently,  other  classes 
of  waveforms  have  been  investigated,  implying  the  possible  addition  of  the 
other  line  segments  shown  in  Figure  6-2. 
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Of  particular  interest  are  the  low  frequency  Crater  and  Related 
Effects  (CARE)  which,  in  Figure  6-2,  dominate  horizontal  response  below 
1/2  cps  and  the  ground  roll  frequency  band  in  which  the  oscillatory  nature 
of  the  groundshock  waveform  interacts  with  shock  isolation  systems, 
causing  amplification  of  peak  relative  displacement.  The  CARE  region 
portion  of  the  spectrum  which  depends  on  the  crater  volume  is  equal  in 
value  for  vertical  and  horizontal  motions,  and  can  disappear  for  non- 
cratering bursts.  All  these  segments  are  related  to  peak  free  field  ground 
shock  predictions  and  the  FAST  code  will  select  random  values  for  these 
variables.  The  proper  way  to  simulate  all  these  regions  of  shock  response 
is  to  input  the  entire  spectrum  for  the  FAST  calculation  and  let  the  com- 
puter determine  the  shock  spectrum  segment  which  dominates  at  the  var- 
ious frequencies  of  interest. 

As  an  example  for  a 600  psi  surface  burst,  consider  the  shock 
spectrum  as  applied  to  the  control  facility's  shock-isolated  floor.  The 
horizontal  natural  frequency  of  the  floor  is  0.  3 cps.  For  5 MT,  this 
frequent  is  in  the  CARE  pseudo- velocity  region  of  the  horizontal  shock 
spectrum.  However,  for  the  low  yields,  at  the  same  pressure,  0.  3 cps 
would  be  in  the  CARE  displacement  region. 

The  d^RE  anc*  VCARE  are  ranc^om  variables,  and  a bivariate 
sample  of  pairs  of  their  values  will  be  generated  during  a FAST  calcula- 
tion. The  bivariate  distribution  will  contain  some  cases  where  the  hor- 
izontal floor  response  is  governed  by  the  leg  a°d  other  cases 

where  the  floor  response  is  governed  by  the  v^^^g  leg  (on  a velocity  leg, 
maximum  relative  displacement  is  inversely  proportional  to  the  frequency, 
i.  e.  , dyej  = V / 2 tt f ).  One  way  to  select  the  proper  value  for  maximum 
relative  displacement  is  to  create  a logical  subnetwork  as  shown  below. 

FLOOR/WALL 

IMPACT 


■ rv^c 
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This  network  will  count  floor  damage  only  if  the  smaller  of 
and  v£ARE^n*h  ifl  large  enough  to  cause  damage.  Note  that  the  concave 
downward  intersection  of  d^ARE  an(^  VCARE  ^nes  results  in  parallel  logic, 
while  a concave  upward  intersection  between  and  the  displacement 

amplification  line  would  require  a series  subnetwork. 

For  height  of  burst  cases,  the  CARE  environments  will  be  small, 
and  the  shock  suspension  system  horizontal  response  will  be  governed  by 
the  air  induced  displacement.  The  situation  where  the  CARE  and  PARE 
environments  are  of  the  same  order  can  also  be  modeled  by  the  FAST  code. 

By  combining  slock  spectra  line  segments  in  series  and  parallel  it 
is  possible  to  model  an  entire  shock  spectrum.  However,  for  the  sample 
problem  system,  soft-mounted  items  respond  only  to  the  lower  end  of  the 
spectrum,  hard-mounted  items  respond  only  at  the  upper  end  of  the  spec- 
trum, and  the  batteries,  which  respond  to  the  mid-frequencies,  have  been 
determined  to  depend  only  on  shock-spectra  velocity.  The  effect  of  scaling 
laws  on  other  components  was  examined  and  it  was  concluded  that  the 
environments  identified  by  the  shaded  numbers  in  Table  6-A  could  be 
dropped  from  further  consideration. 

Derivation  of  covariance  between  environment  systematic  variations 
and  random  variations  was  based  upon  environment  prediction  equations 
using  the  matrix  approach  described  in  Section  2.3.  Table  6-D  shows  the 
starting  point,  which  is  a list  of  the  weapon  effects  environment  equations 
(for  a fixed  range)  in  terms  of  the  analysis  coefficients,  load,  site,  and 
structural  parameters.  The  exponents  of  the  variables  of  the  environment 
equations  are  summarized  in  matrix  form  in  Table  6-E.  In  the  table  the 
environments  are  divided  into  PARE  and  CARE  effects,  with  the  PARE 
effects  further  subdivided  into  free  field  and  structural  effects.  Table  6-F 
shows  the  sources  of  uncertainty  for  the  sample  problem  and  gives  the  K- 
factors  for  both  random  and  systematic  variations.  In  general,  a K-factor 
is  a parameter  such  that  most  random  outcomes  lie  between  K and  1/K 
times  the  best  estimate  value,  M.  K-factors  are  exact  descriptors  of 
variation  for  log-normal  distributions  and  approximate  for  others.  A 
two-sigma  K-factor  corresponds  to  2 a bounds  on  a log-normal  distribu- 
tion; i.  e. , 2-1  /2°7c  of  outcomes  greater  than  KM  and  2-  1 /2%  of  random 
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Table  6-F.  Systematic  and  Random  Variation  of  the  Basic  Parameters 


outcomes  less  than  M/K.  For  this  problem,  2 abounds  are  used.  Hence, 
2 

the  variances,  <rj,  of  the  underlying  sources  of  variation  may  be  computed 
from  the  corresponding  K-factors  as  follows: 


2 ^ 
i “ 4 


(6-4) 


For  this  problem,  it  is  assumed  that  the  random  variation  of  the  analysis 
coefficients  is  small  because  the  analysis  technique  is  based  on  a single 
geology  of  interest  and  the  random  variation  of  the  identified  underlying 
variables  will  account  for  most,  if  not  all,  of  the  random  variance.  There 
fore,  several  of  the  random  K-factors  are  shown  as  unity  in  Table  6-F. 

On  the  other  hand,  the  K-factors  for  systematic  bias  of  the  analysis  coef- 
ficients are  relatively  large,  because  the  trend  lines  of  the  prediction 
equations  are  not  based  on  firm  test  and/or  analysis  data. 

Correlation  coefficients  and  K-factors  for  the  materials,  loads  and 
structural  parameters  are  shown  in  Table  6-G,  including  correlations 
between  M,  c and  p;  correlations  between  p and  I;  correlations  between  G 
and  M;  and  correlations  among  the  structural  load  factors,  where  M is 
confined  modulus,  c is  peak  stress  wave  velocity  and  p is  density  of  the 
soil;  p is  pressure,  I is  impulse,  G£  is  foundation  shear  modulus  of  the 
soils,  and  the  remaining  parameters  in  the  table  are  structural  factors. 


Table  6-G.  Correlation  Among  Basic  Parameters 


Parameters 

P.. 

i 

j 

C1 7 

C18 

0.4 

C12 

C 1 3 

0.  2 

C14 

C16 

0.  3 

M 

c 

0.  8 

M 

P 

0.  3 

c 

P 

0.  3 

I 

P 

0.7 

FR 

FW 

0.  8 

ff 

FW 

0.8 

fr 

ff 

0.7 

Gf 

M 

0.  5 
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The  random  and  systematic  covariance  matrices  of  the  underlying 
variables  (£vv)  can  now  be  constructed  from  the  K-factors  and  correlation 
coefficients.  Then,  /rom  Equation  (2-19),  the  corresponding  environ- 
ment covariance  matrices  can  be  computed.  Table  6-H  gives  2<r  K-factors 
for  the  systematic  and  random  variation  of  the  environments.  The  covari- 
ance matrices  are  contained  in  the  Input  Summary,  printed  out  with  each 
FAST  run.  These  are  in  Appendix  A. 

6.2  TRANSFER  FUNCTIONS  AND  FRAGILITIES 

The  derivation  of  transfer  functions  and  fragilities  is  best  described 
together.  Therefore,  the  discussion  of  both  topics  is  given  in  this  section. 
Two  specific  examples  are  used  to  illustrate  in  detail  the  derivation  of 
transfer  functions  and  fragility  curves,  namely,  for  the  antenna  probe  and 
the  shock  suspension  system  structural  failure.  Derivations  for  other 
components,  while  not  given,  would  follow  a similar  development. 

Sometimes  fragilities  are  defined  in  a normalized  fashion,  by  dividing 
the  unnormalized  fragility  by  the  allowable  load.  In  these  instances,  the 
unnormalized  transfer  function  must  also  be  divided  by  the  allowable  load 
(as  shown  in  figures  of  Section  6.3).  Fragility  and  transfer  function  data 
can  be  displayed  graphically  or  given  in  terms  of  K-factors,  whichever  is 
most  convenient. 

6.2.1  Communications  Antenna  Probe  - Fragility  and  Transfer 
Functloli 

Figure  6-3  shows  a cross-section  of  a communication  antenna  which 
has  been  designed  for  the  hypothetical  control  facility.  The  probe,  of 
width  w and  length  1,  is  cylindrical  as  shown  in  Section  A-A.  The  probe  is 

vulnerable  to  dvnamic  pressure  loading,  which  creates  a moment,  M = 

2 a 
c^  Q wl  /2,  about  the  weak  rectangular  section  shown,  where  is  the 

drag  coefficient  and  Q is  the  dynamic  pressure.  The  exterior  of  the  probe 

is  axisymmetric.  Therefore,  the  applied  moment,  M , does  not  vary  with 

cl 

azimuth.  However,  the  allowable  moment,  M = ir  Z(4>),  is  azimuth- 

P a 

sensitive  because  the  plastic  section  modulus,  Z(4>),  is  a function  of  azi- 
muth, with  the  strongest  section  facing  North,  where  v is  the  allowable 

cl 

stress. 
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Table  6-H.  K- Factors  for  Hostile  Environments 


Environment 

Systematic 

Random 

Number 

Name 

KS 

KR 

1 

d 

z 

2.5 

2.4 

2 

V 

z 

1.9 

2.  3 

3 

a 

z 

2.0 

1.8 

4 

V 

z 

2.0 

2.  3 

5 

20  d 2/v 
z z 

5.0 

2.6 

6 

dH 

2.9 

2.  5 

7 

VH 

2.0 

1.  1 

8 

VH 

2.  1 

1.  1 

9 

1-5vH2/dH 

5.7 

2.  3 

10 

"•  drel 

3.3 

5.2 

11 

Q 

1.6 

1.  3 

12 

aST  ( A / C ) 

1.6 

1.2 

13 

aSTR(SSS) 

1.6 

1.2 

14 

Equivalent  Static  Roof 
Load 

1.4 

1.2 

15 

Equivalent  Static  Wall 
Load 

1.6 

1.3 

16 

Equivalent  Static  Footing 
Load 

1.6 

1.3 

17 

dCARE 

3.0 

1.3 

18 

VCARE 

2.9 

1.8 

19 

Debris  Depth 

2.8 

1.  5 

WEAK  SECTION 


M = APPLIED  MOMENT  = c , Q — 
a d 2 


SIDE  VIEW 


TOP  VIEW  - SECTION  AA 
N 


WEAK  SECTION 

Mp  = PLASTIC  MOMENT  = °gZ(9) 

Z (0)  = 1/4  bh2  (1-1/3t72)/co$0  0 19*37° 


WHERE  n = v/(£cot0)2  + 3 - £cct  9 


Z (4>)  1/4  hb2  ( 1 - y ) / cos  ( 90°  - 4> ) 37°  < <t>  < 90° 


(9)  = 


Figure  6-3.  Antenna  Structure  Description  (Plain  View  and  Elevation) 
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The  critical  response  parameter  is  the  ratio  of  applied  moment  to 
allowable  moment, 


cdwi 


2Z{4>) 


^NORTH 

Z(4>) 


cdwl 


2ZNORTH  'a 


(6-5) 


where  Z^qj^^j  is  the  plastic  section  modulus  versus  north  attack,  and 
ZNORTH^Z^  t^ie  az*mutk  sensitivity  which  is  incorporated  into  the 
transfer  function.  Note  that  the  dashed  lines,  which  are  input  to  FAST, 
accurately  represent  the  continuous  function  both  in  the  mean  and  in  the 
variance  of  the  transfer  function  from  0°  to  360°. 


Azimuth  sensitivity  also  can  be  caused  by  non-uniform  rattlespace. 
Obviously,  non-uniform  rattlespace  is  not  good,  but  sometimes  it  cannot 
be  avoided.  Figure  6-4  shows  a plan  view  of  the  rattlespace  and  the  azi- 
muth sensitivity  which  could  result  for  pure  translational  motion.  In 
practice,  pure  translational  motions  are  unlikely  and  cusps  on  azimuth 
sensitivity  curves  will  be  rounded  by  random  variation. 


6.2.2  S Structural  Failure  Fragility  and  Transfer  Function 


The  following  discussion  considers  only  the  vertical  input  motions 
and  responses,  but  similar  concepts  can  be  applied  to  other  directions. 
The  shock  suspension  system  (SSS)  will  fail  structurally  when  the  floor 
motion  exceeds  the  length  of  the  allowable  excursion  for  the  rod  which 
attaches  the  floor  to  the  isolators.  The  isolators  are  carefully  calibrated 
to  fail  all  at  once  or  not  at  all;  hence,  one  displacement  fragility  curve 
applies  to  all  isolators. 


The  vertical  natural  frequency  of  the  S is  2 cps.  In  some  situations 
2 cps  corresponds  to  the  constant  displacement  portion  of  the  shock  spec- 
trum (see  Figure  6-5).  In  other  instances,  2 cps  intersects  the  constant 
velocity  portion,  and  sometimes  the  displacement  amplification  region  of 
the  spectrum  covers  2 cps.  Transfer  function  uncertainties  are  used  to 
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allow  for  the  limitations  of  the  data  base  and  for  the  straight-line 
simplifications  of  this  possible  cause  of  excessive  excursion.  The  S' 
shock  spectrum  value  may  be  any  of  the  following: 


Portion  of 
Shock  Spectrum 

Peak 

Free  F.eld 
Value 

Constant 

d 

Displacement 

Displacement 

d 

Amplification 

Constant 

V 

Velocity 

Shock  Spectrum  Value 


In  Terms  of  In  Terms  of 
Displacement  Velocity 


Applicable 
T ransfer 
Function 


V/2irf 


20f d /v 


D x 2irf 


The  logic  of  S failure  is  series-parallel.  The  series  portion 

3 ' • 

implies  that  the  S will  fail  if  the  larger  of  the  displacement  amplification 


and  the  peak  free  field  vertical  displacement  exceeds  the  available  stroke 

w 

length.  The  parallel  portion  implies  that  the  subsystem  will  survive  if 
either  a)  the  relative  displacement  at  2 cps  on  the  shock  spectral  velocity 


line  is  less  than  the  stroke  length  or  b)  d^  is  less  than  the  stroke  length 
of  the  S^. 


Figure  6-5  shows  the  effects  of  random  variation  on  the  shock  spec- 


tra. The  elliptical  contours  show  how  one  of  the  shock  spectral  corners 

2 

might  vary.  Also,  shown  are  the  marginal  densities  of  D,  V and  20d  /v 


which  are  all  computed  from  only  2 free  field  variables.  Systematic 
variations  among  these  variables  are  not  linearly  dependent  in  log  space 
due  to  variation  of  analysis  constant  (shown  in  Table  6-F),  which  is  the 
systematic  variation  K factor  for  the  fourth  analysis  coefficient. 


The  cut-out  in  Figure  6-5  shows  the  random  bivariate  distribution  of 
fj*  versus  f^,  which  are  the  random  parts  of  the  shock  spectra  break  fre- 
quencies. Note  that  the  fourth  analysis  constant  in  Table  6-F  has  a random 
K factor,  = 1.  This  implies  that  f *s  constant  and  the 
lower  corner  frequency  is  linearly  dependent  on  the  upper  corner 


■ 
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frequency  (in  log  space).  This  will  result  in  a singular  covariance  matrix 
for  random  environments.  But  this  does  not  matter  since  FAST  can 
decompose  and  use  singular  covariance  matrices  for  random  number  gen- 
eration in  the  simulation. 


6.  3 SYSTEM  NETWORK 


The  preliminary  baseline  form  for  the  system  network  is  shown  in 
Figure  6-o.  This  network  is  the  result  of  a good  deal  of  pruning  of  larger 
networks  based  upon  knowledge  of  environment  distributions  and  component 
criticality.  Each  of  the  subsystem  networks,  together  with  fragility  and 
transfer  function  definitions,  is  shown  in  Figures  6-7  to  6-14 


One  point  not  resolved  in  the  preliminary  system  network,  concerned 
the  air  conditioner.  Two  alternatives  were  postulated  for  the  design  of 
this  subsystem;  one  having  the  fragility  shown  in  Figure  6-14,  the  other 
having  a superhard  design.  Resolution  of  this  decision  was  based  upon  the 
FAST  survivability  statistics. 


6.4  SURVIVABILITY  CALCULATION  RESULTS 


To  illustrate  use  of  the  FAST  methodology  and  code,  assume  a 
requirement  that  the  system  have  a probability  of  survival  of  at  least  0.  8 
with  50%  confidence  when  subjected  to  a surface- burst  nuclear  attack 
involving  a single  600  psi  overpressure  pulse  from  a 5-megaton  weapon, 
The  inputs  derived  in  Section  6.  1 - 6.  3 are  assumed  to  be  for  a baseline 
system,  which  is  the  result  of  an  initial  design  effort.  The  FAST  technique 
is  used  for  evaluating  the  baseline  system  design  and  identifying  how  best 
to  modify  it  to  meet  the  system  requirements. 


FAST  sur\ ivability  statistics  used  to  evaluate  the  baseline  system 
show  that  it  fails  to  meet  the  survivability  requirement.  The  calculation 
also  provides  information  on  how  the  baseline  system  design  could  be 
improved  because  critical  subsystems  and  components  are  identified.  A 
sensitivity  analysis  can  be  done,  employing  a steepest  ascent  approach  to 
identify  the  most  cost  effective  way  to  strengthen  the  system  design  to 
meet  the  system  requirement.  The  steepest-ascent  technique  is  not  an 
inherent  part  of  the  FAST  methodology,  but  is  indicative  of  how  the  FAST 
results  can  be  employed  to  provide  data  for  management  decisions. 


SUBSYSTEM  SUBSYSTEM  SUBSYSTEM  SUBSYSTEM  SUBSYSTEM  SUBSYSTEM  SUBSYSTEM 


Figure  6-8.  Shock  Suspension  System 


Battery  Supports 


tt  :ry  Elements 


Figure  6-11.  Communication  Antenna  System 


Figure  6-14.  Air  Conditioner 


The  sensitivity  analysis  leads  to  a configuration  satisfying  the 
system  requirement  which  is  then  used  to  illustrate  the  effects  of  weapon 
yield  and  height- of- burst  perturbations  on  the  survivability  statistics. 


6.4.1  Evaluation  of  the  Baseline  System 

Evaluation  of  the  baseline  system  was  accomplished  by  means  of  the 
first  computer  run,  which  is  reproduced  in  Appendix  A.  The  inputs  for 
this  run  were  derived  in  Sections  6.  1 through  6.  3. 

An  overview  of  the  results  from  the  first  run  is  given  in  Figure  6-15. 
The  system  probability  of  survival  is  0.05,  much  less  than  the  required 
0.8.  Four  of  the  subsystems  have  a probability  of  essentially  1.00.  Thus, 
to  achieve  the  required  system  hardness  some  combination  of  the  other 
four  subsystems  must  be  hardened.  At  this  point,  it  was  decided  to  use 
the  superhard  version  of  subsystem  8,  the  air  conditioner.  Therefore, 
this  subsystem  was  deleted  from  further  analysis. 

Several  aspects  of  baseline  system  survivability  statistics  are  plotted 
in  Figure  b-16.  The  system  probability  of  survival  systematic  variation 
distribution  (10,  50  and  90%  levels)  is  shown  at  the  upper  left  as  a function 
of  overpressure. 

In  the  upper  right  is  plotted  the  median  value  of  the  probability  of 

survival  of  the  subsystems  as  functions  of  overpressure.  As  can  be  seen, 

3 

the  critical  subsystems  are  the  command  and  control  console  C (SSI), 
shock  suspension  system  (SS2),  and  t'  •»  facility  (SS7). 

The  other  three  diagrams  break  down  the  three  critical  subsystems. 

T le  median  probability  of  survival  of  the  subsystem,  together  with  the 
median  probability  of  survival  of  the  components,  are  shown  as  functions 
of  overpressure. 

Evaluation  of  the  subsystem  survivability  statistics  in  Figure  6-16 
identified  three  areas  where  hardening  is  required.  The  first  is  com- 
ponent C101  in  Subsystem  1,  which  causes  almost  all  of  the  failure  of  this 
subsystem.  The  secoi.J  area  is  Subsystem  2,  the  Shock  Suspension  Sys- 
tem, where  it  is  natural  to  harden  the  three  fragilities  all  in  the  same 
proportion  by  redesigning  the  subsystem.  The  third  area  is  the  facility, 
Subsystem  7.  Here  also,  it  is  natural  to  improve  all  fragilities  in 
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proportion  because  this  is  accomplished  by  using  thicker  concrete.  Of 
course,  in  the  case  of  a real  system,  other  factors  should  also  be  con- 
sidered before  makir.g  a system  change,  such  as  how  well  the  inputs  for  the 
critical  subsystems  are  known,  and  whether  or  not  further  study  in  these 
areas  might  reduce  the  apparent  criticality  of  these  subsystems. 

At  this  point,  three  areas  have  been  identified  requiring  hardening 
to  achieve  the  system  survivability  requirements.  In  the  next  section, 
these  are  treated  as  three  variables  and  a steepest-ascent  approach  is 
employed  to  accomplish  system  hardening  trade  analysis. 

6.4.2  System  Hardening  Trade  Analysis 

Once  one  has  the  three  variables  with  the  predominate  effect  on  sys- 
tem hardness,  the  next  step  is  to  perform  trade  analysis  and  find  a system 
design  satisfying  the  survivability  requirement.  We  shall  neglect  such 
possibly  important  considerations  as  system  costs,  accuracy  of  the  data 
base,  and  confidence  ir  achieving  technical  objectives.  These  factors 
would  have  to  be  considered  in  an  actual  situation  but  are  not  necessary  to 
illustrate  the  concept.  We  will  employ  an  optimization  technique  known  as 
the  method  of  steepest  ascent  to  illustrate  how  such  an  optimization  tech- 
nique can  be  used  with  FAST. 

The  method  of  steepest  ascent  uses  a statistically-designed  experi- 
ment to  find  how  the  independent  variables  (three  hardness  variables  in  our 
case)  affect  the  dependent  variable  (system  survivability  in  our  case). 

This  information  is  then  used  to  compute  the  way  to  simultaneously  change 
the  independent  variables  (called  the  direction  of  steepest  ascent)  to  max- 
imize the  rate  of  increase  of  the  dependent  variable.  System  configura- 
tions along  this  line  are  evaluated  to  explore  the  survivability  response. 
Typically,  several  cycles  of  experimental  designs  and  steepest-ascent 
evaluations  are  employed  to  account  for  non-linearities  and  other  con- 
tingencies. We  shall  limit  our  efforts  to  the  first  cycle,  settling  for  the 
first  system  identified  which  satisfies  the  survivability  requirement. 

The  experimental  design  is  shown  in  Table  6-1.  Two  hardness  levels 
are  used  for  each  of  the  variables,  the  baseline  hardness  and  a higher  value 
by  8 percent  for  the  command  and  communication  console,  15  for  the  sus- 
pension system,  and  8 for  the  facility,  respectively.  Since  there  are 
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two  levels  for  each  of  the  three  variables,  the  total  number  of 
3 

combinations  is  2 =8,  so  this  is  the  number  of  systems  defined  and 

evaluated  by  the  FAST  code  in  the  experimental  design. 

The  input  section  of  the  FAST  printout  for  this  computer  r m is 
reproduced  as  Run  2 in  Appendix  A.  The  probabilities  of  survival  for  the 
eight  system  configurations  are  shown  in  Column  5 of  Table  6-1,  together 
with  the  statistical  analysis  of  the  data. 

From  the  data  in  Table  6-1,  the  direction  5f  steepest  ascent  can  be 
established.  This  was  done,  and  approximately  equally  spaced  points  were 
selected  along  the  line  a's  shown  in  Table  6-J  which  defines  5 system  con- 
figurations along  the  line  of  steepest  ascent.  Along  this  line,  if  fragility 
C101  is  increased  25%,  then  the  shock  suspension  system  fragilities  are 
increased  55%  and  the  facility  fragilities  are  increased  30%. 

The  final  column  of  Table  6-J  gives  probabilities  of  survival  for  the 
five  system  configurations.  As  can  be  seen,  systems  SYS15  and  SYS16 
exceed  the  required  80%  probability  of  survival.  In  the  evaluation  of  a real 
system,  many  more  tradeoff  calculations  would  be  done,  to  fine-tune  the 
design.  But  for  the  purpose  of  this  sample  problem,  system  SYS16  will 
be  selected  as  the  final  configuration.  This  system  configuration  is  used 
in  the  calculations  of  the  next  section. 

6.  4.  3 Yield  and  HOB  Variations 

The  purpose  of  this  section  is  to  demonstrate  the  use  of  FAST  code 
for  calculations  involving  changes  of  weapon  yield  and  height  of  burst 
(HOB).  To  do  this,  the  "optimized"  system  obtained  in  the  last  section 
was  used.  Four  combinations  of  yield  and  HOB  were  employed  as  follows: 


HOB, 


(W/l  MT) 


YIELD,  MT 

0.  5 

5.  0 

X 

X 

X 

X 

Table  6-J.  Results  Along  "Line  of  Steepest  Ascent 


The  input  section  of  the  printout  for  this  computer  run  is  reproduced  in 
Appendix  A.  The  survivability  statistics  for  the  four  cases  are  plotted  in 
Figure  6- 17. 

At  a scaled  HOB  of  1500  ft/(W/l  MT)*^,  air-induced  effects  on  the 
sample  system  are  virtually  the  same  as  for  a surface  burst.  However, 
the  pressure-range  relation  in  FAST  has  made  the  correction  for  the 
ranges  at  which  the  HOB  pressures  are  obtained.  Cratering  effects,  on  the 

other  hand,  are  drastically  reduced.  Crater  volume  scaling  E = E exp 

1/3  ° 

(-X.HOB/W  ) accounts  for  almost  all  of  the  reduction,  while  the  slightly 

larger  scaled  range  for  the  specific  pressure  levels  at  a given  HOB 

decreases  the  crater  induced  effects  by  a negligible  amount.  Similar 

scaling  can  be  applied  for  ranges  where  HOB/R  s 0.4.  However,  this 

scaling  procedure  breaks  down  as  HOB  approaches  the  optimal  HOB  for 

a given  overpressure.  Other  scaling  procedures  based  on  site  response 

to  changing  pressure-time  histories  are  known  to  the  weapons  effects 

community.  These  scaling  laws  can  be  modeled  in  FAST  by  modifying 

environment  and  transfer  function  inputs, 
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